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ABSTRACT 


For the Skylab S-193 radar altimeter, this report presents data 
processing flow charts and identification of calibration requirements 
and problem areas for all presently defined S-193 altimeter experiments , 
and an analysis and simulation of the relationship between one particular 
S-193 measurement ("mean-of-the-peaks") and the parameter of interest 
(”peak-of-the-mean") for determining the sea surface scattering cross- 
section, 0 *^. For the GEOS-C radar altimeter, this report presents results 
of system analyses pertaining to signal-to-noise ratio, pulse compression 
threshold behavior, altimeter measurement variance characteristics, desir- 
ability of on-board averaging, tracker bandwidth considerations, and 
statistical character of the altimeter data in relation to harmonic analysis 
properties of the geodetic signal. 
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CHAPTER 1 


INTRODUCTION AND SUMMARY OF RESULTS 
1*0 SCOPE OF THE REPORT 

This is the final engineering report on Tasks B and E of NASA Contract 
No. NAS6-2135. Task B consisted of a study of data processing and calib- 
ration procedures for the Skylab S-193 altimeter, and Task E consisted of 
system analyses and related factors for the GEOS-C satellite program. Task 
D, which is still in progress and is to be reported separately, comprises 
a study of oceanographic current sensing using remote sensing methods. 

Chapter 2 of this report summarizes our work on Task B, the S-193 task, 
examining the data processing required to calibrate and/or correct the data 
which the S-193 radar altimeter will obtain. Special problem areas such as 
signal- to-noise ratio determination and RF system losses are addressed in 
detail. A data flow and processing diagram which suggests data reduction/ 
analysis responsibilities is also presented. Finally, the problem of relating 
the altimeter AGC measurement to the peak-of-the-niean return waveform is 
investigated in order to obtain radar cross-section (a^) data from the S-193 
altimeter. 

Chapter 3 reports the Task E (GEOS-C) work and gives results of inves- 
tigations of several system engineering aspects of the altimeter. Recommend- 
ed calibration data requirements to support waveform experiments are also 
reported. The system analyses pertain to potential problem areas of signal- 
to-noise ratio (SNR), pulse compression threshold behavior, altimeter measure- 
ment variance characteristics, desirability of on-board averaging, tracker 
bandwidth considerations as influenced by waveform experiments , and the 
statistical character of the altimeter data in relation to harmonic analysis 
properties of the geodetic signal, 

1,1 Conclusions and Recommendations 

The following paragraphs summarize the principal findings of this study 
and the remainder of this report provides the detailed work from which these 
conclusions and recommendations are derived. 

TASK B , the S-193 task (see Chapter 2) : 

1 . Data obtained by the S-193 radar altimeter do not permit a direct 
determination of receiver signal-to-noise ratio. Analytical results presented 
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in this report provide a means for determining the signal-to-noise ratio , 
but they must be substantiated by measurements on the S-193 breadboard 
hardware. 

2. Determination of o° from the S-193 altimeter AGC data can only be 
approximately accomplished for off-nadir angles of 0, 0.5 and 1.5 degrees. 

To overcome this problem, we recommend that an r-f actor study be conducted 
using the S-193 breadboard hardware in conjunction with pulse shaping 
circuits which approximate computed waveshapes for 0, 0.5, 1.5, 3, 8 and 
15.6 degrees off-nadir. 

3. Time sidelobes due to the pulse expansion/ compression network 
prevent the accurate measurement of sample-and-hold (S&H) gate offset levels 
in Mode V of the S-193 radar altimeter. Measurements on the breadboard 
hardware should be made to ensure that S&H level offset is stable from mode 
to mode. 

4. There is a definite need for measurements of instantaneous and 
average tracker jitter as a function of signal-to-noise ratio. These data 
are required for post-flight S&H gate position reconstruction and evaluation 
of the altimeter under radial acceleration conditions. 

5. Additional RF-loss measurements must be made on the breadboard 
hardware to eliminate unknown loss factors which directly affect the measure- 
ment of a . 

TASK E , the GEOS-C task (see Chapter 3) : 

1. We recommend the incorporation of a lower (^.5 Hz) closed-loop 
altitude tracker bandwidth option in addition to the presently planned 4.0 Hz 
bandwidth. The availability of this lower bandwidth option would lessen the 
need for derivation of S&H positional jitter in waveform reconstruction, 
provide more desirable slew characteristics for the final acquisition step, 
and provide an on-board filtering capability more nearly in consonance with 
the spatial filter effect. More study and analysis is needed to quantitatively 
characterize this last effect, especially in regard to data processing methods 
for altitude and slope extraction. 

2. Regarding system configuration, we feel that the on-board averaging 
processes (the signal conditioners following the altitude tracker) do not favor- 
ably affect statistical properties of the altimeter data and are generally 
superfluous. These characteristics may in fact be obtained by appropriate 
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changes in tracker bandwidths. Also signal-to-noise considerations are felt 
to argue against use of a 40 dB time sidelobe filter with its attendant loss 
properties. A 20 dB filter is considered to be a better choice. Also, SNR 
levels are such that time sidelobes are not expected to be a prime consideration. 

3. We reconmiend that system level tests of the type described be 
conducted to determine if latent problem areas exist, such as the threshold 
effect discussed. We also emphasize the need for statistical data on the 
altitude quantization problem and on tracker performance at SNR values near 
0 dB. 

4. The acquisition of more performance and qualification test data is 
strongly recommended. The Skylab S-193 altimeter program has amply demonstrated 
the need for such data in support of data processing and technology evaluation. 
We also recommend that the GEOS-C Altimeter Experiment Proposals be reviewed 
before completion of the hardware test phase to determine the possible need for 
added test activities, even if the satellite program time schedule dictates 
that such data be obtained on back-up hardware. 
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CHAPTER 2 


TASK B: S-193 ALTIMETER CALIBRATION LOGIC 

2.0 INTRODUCTION 

As a consequence of the spacecraft hardware constraints and the rather 
sophisticated experimental objectives of the S-193 radar altimeter, little if 
any of the raw data from the altimeter can be used without extensive calibra- 
tion. Except for the Nadir Align mode, each functional mode of the altimeter 
contains calibration steps in its operating sequence. It would of course be 
desirable to obtain all of the calibration data during in-flight operation of 
the altimeter; however, this is impractical. Therefore, it will be necessary 
to augment the in-flight data with preflight measurements. In addition, some 
calibrations cannot be accomplished even during preflight, and we must rely 
on simulation and/or analytical results. A very general block diagram of how 
these various calibration data will be merged to yield information for the 
sensor evaluation and EREP programs is shown in Figure 2.1. 

Except for the requirement of proper operation of Mode VI (Nadir Align) , 
al-l of the operating modes of the S-193 radar altimeter are functionally 
independent. Also, the data obtained in one mode may not be completely cor- 
related with the data from another mode due to the motion of the spacecraft 
and changes in the scattering surface. For this reason, it appears logical to 
develop raw data calibration procedures on a mode-by-mode basis . Basic data 
calibration block diagrams are developed in Section 2.1 of this chapter. 

Once the raw data are calibration corrected, they are available for use by 
S-193 radar altimeter sensor technology and EREP experimenters. The goals of 
these experimenters, however, dictate additional data processing which is 
beyond mere calibration. In order for NASA to know the type and scope of this 
additional data processing. Section 2.2 presents a general discussion of this 
topic. In particular, a block diagram' is presented which outlines the antici- 
pated areas of responsibility and the type of processing anticipated for the 
EREP categories. Also shown in the diagram are the expected outputs from the 
various experiments. Additionally, Section 2.2 contains a short discussion of 
the type of data reduction which will be necessary in order to obtain a "quick- 
look" estimate of o° from the Skylab data. 
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2.1 Altimeter Data Calibration by Mode 

2.1.1 Mode I; Pulse shape .— Throughout this section we will refer to the 
mean location of the returned signal relative to the noise, ramp and plateau 
gates and also to the locations of the Sample and Hold (S&H) gates in the various 
operating modes. To avoid possible confusion, Figures 2,2 and 2.3 have been in- 
cluded to show the positions of the various gates relative to the mean return 
for all of the operating modes and DAS (Data Acquisition Step) or CDS (Calibra- 
tion Data Step) sub-modes. Figure 2.2 is applicable to Modes I, II, and V, while 
Figure 2.3 covers Mode III. In Figure 2.3, the time difference between the 
first and second returns is not to scale in order to get all of the 2nd returns 
on the graph. Also, the number "4" appearing in the S&H gate positions in 
DAS-2 and 3 denotes the presence of four S&H gates; the relative spacing of 
the gates is the same as the S&H gates which sample the first return. 

Figure 2.4 is the data processing block diagram for obtaining calibration 
data from Mode I, CDS-1, 2, and 3. The purpose of. this processing is to obtain 
corrected data which are not available from the DAS sub-modes. In particular, 
the corrected outputs will bei the peak transmitter power, S&H gate timing, 

S&H gate level offset and transmitted waveforms as influenced by the transmitter, 
receiver and S&H gates. Most of the processing indicated is relatively straight- 
forward; however, there are some special points of caution to be noted. Level 
offset of the eight S&H gates should be determined from the first eight S&H gate 
positions in CDS-1 since the gates are sampling zero-mean noise in this sub-mode. 
It is also recommended that a level offset be determined from the last eight 
S&H positions in CDS-3 (see Figure 2.2), and these results should also be the 
same as with CDS-1, Histograms should be plotted for these S&H gates to 
determine if each S&H gate is operating properly.. Any drastic difference in 
the histograms recorded by each of the eight S&H gates (apart from changes in 
the mean) will most probably indicate a malfunction in the gate for certain 
classes of malfunctions. In this way we can isolate improperly operating 
gates and ignore the data recorded by the gate. L* and L* are losses in the 
transmit and receive paths, respectively, which are not included in the path 
loss measurement made in the CDS sub-modes. A more detailed discussion of these 
losses is presented in Section 2.1.5. 

Figure 2.5 indicates the processing which will be required for all of the 
DAS sub-modes of Mode I. As shown by the stars, there are a number of data 
processing functions which are either yet to be determined or are problem 
areas. Included in this category are the following outstanding problems. The 
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Figure 2.3. S&H gate positioning for Mode 3 based on information available as of November 1972. 




Figure 2,4. Data processing diagram for Mode I, CDS-1, 2, and 3. 
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Figure 2.5. Data Processing for Mode I, DAS-1, 2, and 3. 











ratio of the peak-of-the-mean return waveform (which is required to determine 
cr*^ and Signal-to-Noise Ratio (SNR)) to the mean-of-the-peaks of the return 
waveform (which is measured by the AGC) is defined as the r~f actor . Determina- 
tion of the r-f actor is a major problem and is discussed in Section 2.1.5 and in 
Appendix 2B. Compensation of the tracker jitter in order to align the S&H 
gates appears to be a very complex problem because of the different signal 
points involved and the marginal rate at which the tracker error is recorded 
(eight times per second) . Details of a tracker simulation scheme are given 
in Appendix 2C. It is also not obvious how an estimate of the SNR can be made 
from the altimeter data. Originally, it was thought that the integrated noise 
gate voltage would be a valid measure of the post-detection noise power. 

However, since the input to the noise gate integrator is AC-coupled, this 
clamps the integrated noise gate voltage to zero. Thus, determination of the 
SNR is not straightforward and we discuss it further in Section 2.1,5. It 
should be noted (see Figure 2,5) that we are suggesting that the fluctuation 
statistics be computed in the plateau gate region of the mean return. The 
reason for selecting this region is that the mean return waveform is least 
rapidly varying in the plateau region (see Appendix 2A for mean waveform 
computations) and the S&H gate position uncertainty relative to the mean return 
waveform will least affect the fluctuation statistics here. Also indicated on 
the DAS diagram for Mode I (Figure 2.5) is the fact that the S&H offset data 
and the peak transmit power are derived from CDS data. 

2.1.2 Mode II: Cross-section (a°) .~^The processing of the CDS data from 

Mode II is exactly the same as for the CDS data from Mode I (see Figure 2.4). 

To understand the processing that must be accomplished for the DAS sub-modes 
of Mode II, it is first necessary to understand how is to be extracted from 
the data. In Appendix 2A, it is shown that may be computed from a knowledge 
of the peak-of-the-mean return which depends on the antenna gain, pattern and 
pointing angle, altitude, system and propagation* losses and transmitted wave- 
form and power. Thus, if the altimeter measured the peak-of-the-mean return 
and we knew these other parameters, we could compute from the measured 
data. All of the above parameters can be measured except the peak-of-the-mean 
return. The actual S-193 hardware measures the mean-of-the-peaks instead of the 
peak-of-the-means , and thus we must obtain a relation between these two 
quantities, i.e., the r-f actor. This particular problem is discussed more 
thoroughly in Section 2.1.5 and Appendix 2B. For the present purposes, 

*Propagation losses over and above free-space losses. 
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assume that this problem can be resolved and that we know, for eas-h pointing 
angle, the relation between the peak-of-the~mean and the mean-o.-the-peaks 
or the r-factor. The data processing diagram in Figure 2.6 illustrates the 
basic operations which must be accomplished to obtain a from the data. In 
DAS-2 through 6, the tracker will not be operating, thus there is no way to 
compare the computed and measured waveforms. Other than obtaining a relation 
between the peak-of-the-raean return waveform and the mean-of-tbe-peaks (the 
r-factor) , data processing for this mode appears to be relatively straight- 
forward. 

2.1,3 Mode III: Time correlation .— The processing for the CDS sub-mode 

data from Mode III is shown in Figure 2.7 and is essentially the same as for 

Mode I except that the S&H gate and pulse-to-pulse spacing must be accurately 
known. There is a problem with determining the level offset of tlie S&H gates. 
Figure 2.3 shows that in CDS-1 and A, only the first four or five S&H gates 

will sample the noise prior to the CDS pulse. Thus, we will not be able to 

determine a level offset for the last three or four S£^H gates. Hopefully, the 
level offset will be stable enough' that offset data from Modes I or II can be 
used for l*lode III. 

The general scheme of data processing for the DAS sub-modes of Mode III 
is shown in Figure 2.8. This mode was designed to provide the opportunity 
to study pulse-to-pulse correlation as a function of pulse spacing, SNR and 
pulse length. The transmitter, receiver and scattering surfaces are the 
primary contributors to pulse-to-pulse correlation effects. Irora the data 
available from Mode III, it is not possible to isolate the effects of each of 
the above factors. From power spectral density computations given in ref, 2.5, 
we estimate the receiver decorrelation time to be on the order of 15 to 20 ns 
and 150 to 200 ns for the 10 ns and 100 ns pulse widths, respectively. 
Unfortunately, there does not appear to be any way to verify these figures 
from the data available from Mode III; however, since the minimum pulse 
separation in the dual pulse configuration is one microsecond (ps) , we anti- 
cipate that any observable pulse-to-pulse correlation will be a result of 
scattering surface effects. 

In all of the dual-pulse sub-modes of Mode III, four of the S&H gates 
sample the first pulse, while the other four S&H gates sample the second 
pulse. After one frame of data (1.04 sec, or about 250 pulse-pairs), the 
S&H gates are shifted to sample another portion of the return pulses. After 
sequencing through the appropriate number of S&H gate positions (four for DAS-2, 
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Figure 2. 6, Data processing diagram for Mode II, DAS 1-7: note that DAS 2-6 will not 
provide S&H data. 
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Figure 2.7. Data processing diagram for Mode III, CDS-1, 2, 3, 4, and 5. 
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eight for DAS-3) for each sub-mode, the pulse-pair spacing is changed. Thus, 
for each S&H gate position and each pulse-pair spacing, we will have about 
250 returns to work with. 

A typical sequence of return pulses is shown in Figure 2.8a. The S&H 
gates are represented by the vertical arrows and are shown to be in the first 
position (see Figure 2.3 for other typical positions). The pulse-pair spacing 
is T^ and the pulse-pair period is T (4 ms) . 

The following notation and discussion are restricted to a single position; 
the measurement process described is to be repeated for each of the other 
positions as well. The (local) time origin is taken as zero at the first S&H 
gate in the first sample set after the position is set. Then the voltages 
recorded by the S&H gates for the first pulse and the first sample set are 
(for T = S&H gate spacing) , 

Vj^(O) , V2 (t), V2(2t), v^(3t), 

while the voltages recorded for the S&H gates in the second pulse of the first 
sample set are 

v^(y, y{2x + t^). + t^) . 

The voltages recorded by the S&H gates for the first and second pulses of the 
ith sample set are, respectively, 


v^((i - 1)T), V 2 (x + (i - 1)T), V3(2t + (i - 1)T), v^(3t + (i - 1)T) , 

and 

v-(T + (i-DT), v.(t + + (i-l)T), v„(2t + T^ + (i-l)T) , v,(3t + T^ + (i-l)T) . 

±p ‘t p 


Because of the paucity of S&H values for each pulse, it is not advisable to do 
a conventional cross-correlation of the first and second pulses. What is really 


required is the correlation between corresponding S&H gates in the first and 
second pulses of the pulse-pair. That is, we need to know the correlation 


between v^(0) and v^(T^) , v^(t) and ^^(t + T^) , etc. 

If we define p(Jl, T^) as the cross-correlation coefficient* of the £th S&H 
gate for a pulse-pair spacing of T^, we have 


*p = 0 implies independence, while p = 1 implies complete dependence. 
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Figure 2.8a. Two typical pulse-pair returns with the S&H gates shown in the first position. Note 
that of the eight total S&H gates, four are allocated to the first pulse and four are 
allocated to the second of the pulse-pair. 
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1 
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i=l 
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+ (i-l)T) - V 


.H 
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i p 


(i-l)T)- v^l 


.2 


where is the mean of the £th S&H gate for N measurements, 

N 

= |Ev^((JI-1)t + (i-l)T) 
i=l 

, ~~2 

and is the variance. 


^ = i E + (i-l)T) 

i=l 


N 


- -2 

The second-pulse mean, v , and variance, v , are defined similarly. By 

Xj X/ 

expanding terms inside the summations in the expression for p(£,T^), one may 
obtain the following alternative expression: 

N 


p(A,Tp) = 


I p + Tp + (i-l)T)| 


2 -2 

V V 

a £ 


^ )E v,((^-l)t + (i-DT)j js v,(a-l) 


T + Tp + (i-l)T)| 


2 ^2 


This expression for p(£,T^), for £ = 1, 2, 3, or 4, is somewhat more computa- 
tionally convenient to use than the earlier p(£,T^). 

In these expressions, N is the total number of pulse-pairs sampled ('^250) 
at a given position and pulse-pair spacing. For a given pulse-pair spacing, 
p(£,T ) should be independent of £; that is, each point on the first pulse and 

K V ^ P 

each point on the second pulse which are separated by should yield the same 

cross-correlation coefficient, independent of jitter. The results of this 

/ 

measurement and calculation process for a given position of the S&H gates will 
allow us to determine p(£,T^) as a function of T^ and consequently to determine 
how the scattering surface decorrelates. As previously noted, the above 
description has been for a single position, and carrying out the same procedure 
for the other positions will yield additional estimates of p. 
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2.1.4 Mode V: Pulse Compression . —As shown in the CDS data processing 

diagram for Mode V (Figure 2.9), the primary emphasis in this m )de is on 
obtaining data on the pulse shape after the expansion/compression process. 

The reason for this concern is that we do not want to confuse hardware problems 
such as timc-sidelobes with sea-state effects. In this study, it would have 
been desirable to compute mean return waveforms for the 10 ns pulse compression 
case as was done for the 100 ns pulse length case (see Appendix 2A) , since we 
would then have some estimate of the effect of the time sidelobes on the mean 
return waveform. However, time was not available to do this. If relatively 
high time-sidelobes do appear in the pulse as recorded by the CDS -2 sub-mode , 
it is not entirely clear at this time how best to correct the data from the 
DAS-2 sub-mode to eliminate this effect. The data currently available on the 
pulse compression network do not appear to be extensive enough to obtain the 
transfer function of the system. Therefore, it does not seem feasible to 
consider a deconvolution process for correcting the data. 

The data processing for the DAS sub-modes of Mode V is diagrammed in 
Figure 2.10. As with some of the other modes, the problems of SNR and tracker 
variance determination are present. Obtaining the proper level offset compen- 
sation of the S&H gates for all of the DAS sub-modes may be difficult because 
time sidelobes may fall -into -the- -region- where we would normally expect just 
noise (see Figure 2.2). Hopefully, this mode will provide a more favorable 
SNR with which to estimate sea state from the mean return waveform. 

2.1.5 Special Problems to be Resolved .— Throughout the previous sections of 
this chapter we have made numerous references to certain problems which have a 
direct impact on the experimental data that may be obtained from the radar 
altimeter. The purpose of this section is to present detailed discussions of 
these problem areas and to recommend possible solutions. 

SNR Determination : For purposes of design of future sensors, and validity 

of S-193 data, it is most important that the SNR be known in each of the operating 
modes. In order to obtain this ratio, it is necessary to measure the noise power 
in the receiver prior to the arrival of the scattered pulse. Originally, it was 
thought that such information could be derived directly from the S&H gates which 
are positioned in the noise gate. In order to maintain amplifier stability (ref. 
2.1), GF inserted a large blocking capacitor between the output of the video 
amplifier and the S&H gates input. Unfortunately, this capacitor converts the 
noise voltage (at the input to the S&H gates) to a zero-mean process, and the mean 
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Figure 2.9. Data processing diagram for Mode V, CDS-1, 2, and 3. 
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Figure 2.10. Data processing diagram for Mode V, DAS-1, 2, and 3. 










of the S&H gate values is no longer proportional to t’,i noise power. While 
we can construct histograms for the noise voltage recorded by tne S&H gates, 
there is no direct measurement which cap be made to determine the mean noise 
power. 

In order to obtain an approximation to the noise power, we will use the 
results of ref. 2.5. That is, assuming that the noise into the IF filter is 
narrow-band Gaussian with zero mean, we can compute the noise power as recorded 


by the S&H gates. This process will yield a relationship between the post-IF 
noise power and the noise power as recorded by the S&H gates (after the blocking 
capacitor). Through this relationship we can relate the actual noise power as 
measured by the S&H gates to the post-IF noise power. The basic assumption in 
this process is that we know the characteristics of the various receiver com- 
ponents, i.e., the IF and video filters, the detector and the blocking capacitor. 
Thus, we see that accurate measurements of these receiver characterisitcs are 
essential to the qualification of this process. 

The important parts of the receiver are shown in Figure 2,11. It should be 
noted that the IF and video filter bandwidths scale down by a factor of ten when 
switching from 10 to 100 ns pulse length operation, thus the results should apply 
to either case. We assume that the noise into the IF filter is flat over a band- 
width large compared to the IF filter bandwidth and has a spectral amplitude of 

N watts/hertz. We also assume that the noise has a Gaussian density with zero 
o 

mean. The power spectral density (psd) of the noise is taken to be 


N 

o 

$ (to) - 
n 

0 


I to I 

|co| > to^ 


The post— IF psd is given by 

4>IF(iii) = ^‘^(to) |H^p(to)l^ , 


where H (to) is the transfer function of the IF filter. 

IF 

Since is assumed to be much greater than the IF filter bandwidth, 


The noise power out of the IF filter is where 
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The psd at the output of the detector is given by 

= 2a^ I 

where the symbol * denotes convolution. 

The psd at the input to the S&H gates is 

= 2a^ |H^(oi)|^ 

+ 2a^ Njj, |h^(u)|^ |Hcz(“)|^ > 

where H ((a)) is the transfer function of the video filter and is the transfer 

function of the ZC-circuit at the output of the video filter. The noise power 
as recorded by the S&H gates is given by 

00 

= J: — / (J) (qj) (Jq] (2-3) 

SH 2tt f SH'^ ^ 

Zoo 


2 „2 


+ C( N 


IF 
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( 2 - 2 ) 


or 
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SH 


2 2 
2a^N 
o_ 

2tt 


oo 

/I 




‘IF 




dco 


(2-4) 


In obtaining (2-4) we have used the fact that |h^ 2(0) | = 0. Using (2-1), we can 
solve (2-4) for in terms of i.e.. 


00 



Equation (2-5) is the desired relation between the noise power at the S&H gates 
and the post-IF noise power. Its accuracy depends almost entirely upon how well 
we know the various transfer functions. The transfer function about which we 
have no knowledge is H (w) , and we recommend that this be measured on the S-193 

i-iZf 

altimeter breadboard. It should be noted that since the S&H gates are sampling 
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zero mean noise, is just equal to the variance of the voltage recorded by 

b H. 

any one S&H gate in the noise region prior to the arrival of a return pulse. 

RF Hardware Losses : One of the primary purposes of the CDS sub-modes is to 
provide a check on the peak transmitted power. This check is implemented as sho^vni 
in Figure 2.12. A switch between the transmitter and the coupling circulator is 
activated and the transmitted power* is routed through a fixed attenuator and then 
into the front end of the receiver. If the loss in the attenuator is and the 
loss through the switch and waveguide connecting the transmitter (through the atten- 
uator) and the receiver is L^, the peak power into the receiver is P^/(L^L^). Com- 
pensating for the known loss through the attenuator would lead to a measured trans- 
mitter power of Unfortunately, in a DAS configuration the actual loss through 

the waveguide and circulator hardware will be - L* L* (see Figure 2.12). 

In the computations for 0^, we need the ratio where L is the actual 

attenuation between the output of the transmitter and the input to the antenna 

terminals plus the attenuation between the antenna terminals and the receiver 

input, i.e,, L = L* L*. Data from the CDS sub-modes will yield rather 

than P /L* L* as shown above. Thus, there exists a very definite need for a 
t T R 

pre-flight measurement of the quantities L* and L*, for only then can we 
adequately correct the CDS data. It should be noted that any difference between 

and L* L* will be directly reflected in the computation of (see Appendix 2A) . 
a T R _ o 

While this may seem to be a minor point, it could conceivably cause an error in a 

of 1 to 2 dB depending on the ratio LVL* L*. To date, we have seen no GE calib- 
ration or qualification data which provide a measure of any of these required 
losses . 

r-Factor Determination: As we noted in the discussion of Mode II (Section 

2.1.2), before can be extracted from the measured data, an accurate study of 
the relation between the mean-of-the-peaks of the return (AGC voltage) and the 
peak- of-th e-mean return (required to compute a°) must be accomplished. The first 
step in obtaining the relation is to compute the mean return waveforms and, hence, 
the theoretical peak-of-the-mean return. This computation must be accomplished 
because the relation between the two quantities is highly dependent upon the 
amplitude and the shape of the mean return waveform. A computation of the mean 
return waveform is given in Appendix 2A for DAS 1 through 7 of Mode II. 

GE made some preliminary computations of the r-factor (ref. 2.3), but these 
calculations were necessarily approximate since some of the parameters used were 


*0f peak amplitude P^. 
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‘Figure 2.12. S-193 altimeter configuration during measurement of peak transmitted power 

and pulse shape (CDS Sub-mode) . 








not Skylab parameters. We performed an analog/digital computer simulation of 
the altimeter hardware in an attempt to verify GE's results. The results and 
details of the simulation are given in Appendix 2B. Although the simulation 
effort was not complete in all details, it clearly showed that the presently 
available data are questionable , We therefore recommend a much more detailed 
analysis for the determination of the r-factor for all pointing angles of the 
antenna. It is essential that this simulation effort be accomplished before 
the Skylab launch data, for if we do not know the r-factor, we cannot determine 
if Mode II is providing reasonable values of cf°. That is, there will be no way 
we can obtain "quick-look" data on from the Mode II measurements. A block 
diagram of how the r-factor can be determined by simulation is shown in Figure 
2.13, however, hardware measurements are strongly recommended, 

2.2 Altimeter Experimentation Data Flow 

The data processing diagrams presented in Section 2.1 were primarily 
directed toward calibrating the data from the S-193 altimeter. How the S-193 
altimeter data are used and what results are sought are major areas of concern. 

It is also essential that some decision be made on what organization is to be 
respons ible_f or .the. yaripus areas of data processing. 

This section presents a data processing and usage chart which attempts to 
answer some of the above-stated questions. Also presented in this section is 
the processing which must be accomplished in order to obtain a "quick look" value 
for 0° from the data. 

2.2.1 Data Flow and Processing Diagram .— Construction of a data flow block 
diagram for the S-193 altimeter involved establishment of certain ground rules. 
These so-called ground rules were based either on known facts or on our estimate 
of certain situations. Using these criteria, we established the following guide- 
lines for preparing the S-193 altimeter data flow diagram: 

a. data inputs would comprise only the following: 

S-193 altimeter data, 

pre-flight calibration and/or simulation data, 
ground truth data, and 
spacecraft tracking data; 

b. NASA-MSC (Houston) in conjunction with NASA-Wallops Station would: 
verify proper altimeter operation, 

accomplish "quick-look data processing, 
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Figure 2.13. Block diagram of how r-* 
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perforin all engineering units conversion, 
distribute standardized data tapes to experimenters, 
accomplish basic data calibration and correction,* 
collect and refine ground truth data, and 
distribute corrected S~193 and ground truth data to the 
experimenters ; 

c. the S-193 altimeter experimental programs would be: 
sensor technology, 
sea-state estimation (EREP) , 
altimetry, and 
terrain scattering (EREP) . 

In addition, based on available information, we assumed that the outputs from 
the S-193 experimental programs would be as follows; 

a. sensor technology: 

pulse compression evaluation, 
tracker bias determination, 

evaluation of tracker performance versus SNR, 

evaluation of tracker averaging effects, and 

estimation of decorrelation time of the scattering process; 

b. sea-state estimation (EREP): 

0 ° determination, 

_ o 

sea-state estimation from o , 

sea-state estimation from mean waveform smearing, and 
comparison of o° and waveform sea-state estimation techniques; 

c. altimetry; 

geoidal feature extraction, and 

determination of sea-state effects through radar-derived pdf of 
waveheight ; 

d. terrain scattering (EREP) : 

measurement of 0 ° from various terrains , and 

determination of specular scattering effects on the altimeter. 

Figure 2.14 illustrates how and by which experimental program the S-193 
altimeter data will be analyzed. Also sho^ra are the ground truth data require- 
ments which are based on the goals of the EREP programs. The terrain scattering 
and altimeter program were left blank since it is not known how the data are 
to be analyzed. 

*To Include the processing as indicated in the charts in Section 2,1, 
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2.2.2 Quick Look Data Processing ?br \^rtiile most of the data processing 

in Figure 2.14 will be accomplished post-mission, there is one set which must be 
performed soon as the data are available.' In order to determine if the S-193 
altimeter is operating properly and also to plan the next phase of the mission, 
NASA-MSC must have the capability to accomplish "quick-look” processing of the 
data. From our analysis of the situation, Mode II (a°) appears to require the 
least amount of data processing since many of the unknovi?ns can be determined from 
"look-up" tables. In this section we will present a means for obtaining a from 
Mode II data and look-up tables. 

It is shown in Appendix 2A that the peak of the mean return power is given 
by 


<P (T ,^)> = 
^ rss p ^ 




_ ^ H 


L* 

R 


n 

r G X “I 

o 


3 2 

_(4tt) h _ 


F(h, a, Tp) a (C) (2-5) 


where 


T *= Time at which peak of the mean return occurs, 

P 

G = Antenna boresight power gain, 
o 

yv 

= Peak transmitted power, . 

L* L* - RF losses due to the waveguide hardware between the 
T R 

transmitter output and antenna input and the antenna 
output and the receiver input, respectively, 


Lp = Propagation loss over and above free space loss, 

X = RF wavelength, 

h = Altimeter altitude above the surface of the earth 


a = A parameter determined by the width of the pulse 
recorded in Mode II, CDS-1 (~2100) , 

= rms waveheight of the sea surface, 

^ = Pointing angle of the antenna relative to nadir, 

F(h,C ) = An integral which will be provided in a look-up table, 

m p 
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and o°(0 is the desired backscattering cross section at an angle ^ off nadir. 
The digitized AGC voltage may be converted to » ^^peak^ through the use 

of GE-provided look-up tables (ref. 2.4). The peak-of-the-mean return power in 
(2-5) may be computed from 


<^rss<^p’«> = <"rss^^-^>peak> 


( 2 - 6 ) 


where r(0 is known for ^ = 0^, 0.5° and 1.5 from Appendix 2B and as yet is un 

determined for ^ ^ 3, 8 and 15.6°. In Mode II, CDS-1 the AGC measurement will be 

related by GE look— up tables* to (P^/L^L^) where — 119 dB and is a quantity 

yet to be measured (see Section 2.2.5). For the present, we shall assume that 

L* » L* L*, but this assumption must be verified by measurements on the engineering 
a T R 

model. 

The parameters measured in the altimeter and converted by GE look— up tables 


are 


=<"rss<^’^)peak> 


(2-7) 


and 


(M 2 ] 




( 2 - 8 ) 


Substituting (2-8). (2-7) and (2-6) In (2-5) yields the following; 


r(0 Mj^(6) = [M 2 L^l 


g2 xM 

° 3 ' 2 [ ,?,a,T ) o°(0 

(4ir) h ) 


(2-9) 


Solving (2-9) for o°(C) in decibels gives 


[o°(6)]jg - [Mj^(6)]jg 


r(C) 


2 2 
G X 

O n ^ > C > ? » s ^ ) 

( 4 tt) h P 


dB 


*A measurement of receiver temperature is also necessary to insure that 
the correct look-up table is used. 
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The last factor in the above equation can be computed and provided in a look-up 
table, while the first two quantities are the result of measurements (AGC A/D 
converter output) and a look— up table (AGC A/D converter output vs . input power 
and receiver temperature). If we denote the last two terms in the above equation 
by Q, we have 




(2-10) 


The complete processing necessary to obtain "quick-look" values of o is shown in 
block diagram form in Figure 2.15. No requirement for the rras waveheight, has 
been shown since computations (Appendix 2A) have shown that F(h,C ,^,a,T ) is 

iU ^ 

insensitive to C for 0 < C <2 meters. Values for F(h,? niay be deter- 

^m — m — m p 

mined from curves presented in Appendix 2A for an altitude of 235 nautical miles 
and a measured a of 2100. 
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APPENDIX 2A 


MEAN RETURN WAVEFORM COMPUTATION FOR MODE II 

The primary purpose oF Mode II is to provide data on o , the back- 
scattering cross-section of the ocean’s surface per unit area of scatter- 
ing surface. To obtain o° from the data provided by Mode II, it is first 
necessary to compute both the shape and amplitude of the average back- 
scattered pulse. According to the linear scattering theory (ref. 2A1) , 
the average or mean return waveform is a convolution of the following 
functions : 

(1) transmitted pulse waveform, 

(2) flat sea impulse response, 

(3) radar observed wave height distribution, and 

(4) radar receiver impulse response. 

The flat sea impulse response, above, includes the antenna pattern and 
also variations of with the angle of incidence. Since the convolution 
operation is both commutative and distributive, we are at liberty to 
regroup the above convolutions in any desired order. In particular, the 
convolution of (1) and (4) is easily seen to be the waveform as recorded 
by the Sample and Hold (S&H) gates in Mode II, CDS-1. Using the pulse 
shape as measured and calibrated by GE (ref. 2A2, page 6-86) plus the fact 
that both the S&H word count-to-voltage conversion* (ref. 2A2, pages 6-38 
to 6-85) and the detector power-to-voltage conversion are linear, the 
normalized pulse shape was found to be nearly Gaussian and of the form 


P(tp 


P 

t 

L 



e 


(t^ - 100)^ 
a 


(2A1) 


where a = 2100 for the time t^ measured in nanoseconds, P^ is the peak 
transmitted power and L is the system loss. 


*At +25®C and over a word count range of 130 to 70 counts. 
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.For the 100 ns pulse length used in Mode II, it is reasonable to assume 
that the radar observed wave height distribution is Gaussian. Such higher 
order effects as skew and kurtosis may be significant in the case of 10 ns 
pulse scattering, but the 100 ns return is not expected to be sensitive to 
such small pertubations . Therefore for Mode II, the wave height distribution 
p(z), where z is the height above the mean plane (z=0) , is taken to be, 



(2A2) 


where ^ is the mean-square wave height, t - 0 in equation (2A1) corresponds 
to the time at which the leading edge of the transmitted pulse, scattered from 
the surface at a height of h— z beloxi7 the radar, enters the receiver. If the 
variable t denotes the time relative to transmission of the pulse (t=0) , then 
t^ - 0 corresponds to t = 2(h-z)/c. Converting from t^ in- (2A1) to t yields* 


p(t+!) =r 


- i [t - 2(h-z)/c-100]^ 
e a 


(2A3) 


'T\a 


The effect of sea-state roughness upon the transmitted received waveform may 
1)0 -fetmd' from -the-fo flowing convolution: 


P 

ss 


(t) 



P(t + z/c) p(-z) dz 


(2A4) 


Letting x = t -2h/c, i.e., letting x be the incremental two-way ranging time 
relative to the mean height h, and accomplishing the integration in (2A4) yields 



ac 


(2A5) 


*To be consistent with (2A1) , t should be in units of nanoseconds and 
c == 0.3 meters/nanoseconds . 
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We now must determine the flat-sea impulse response and convolve this result 
with (2A5) to find the complete mean return waveform. 

It has been previously shown that the flat-sea impulse response is given 
by ref. 2A3 


r S(t - — ) a°(if) 

° ,o I ^ 7 dA (2A6) 

L (Ait) / r 

scattering 
surface area 

where ; 

’ L = Propagation loss over and above free-space loss, 

= Boresight power gain of the transmitting/receiving antenna, 

X = rf-wavelength , 

2ir 

6(t - — ) = Delta function, 
c 

fC6,oj) = Antenna voltage pattern normalized to a maximum amplitude of unity, 

r = Distance from the antenna to the incremental area dA on the 
scattering surface, and 

= Backscattering cross-section of the sea surface per unit area 
of scattering surface. 

The angle between the boresight axis of the antenna and the nadir axis is 
while T is the aximuth angle on the scattering surface of the antenna 
boresight axis. 9 is the angle between the boresight axis and a line extending 
from the antenna to the elemental surface area dA, and m is the projection of 
the pattern aximuth angle on the scattering surface. The other various geo- 
metrical parameters are defined in Figure 2A1. The antenna pattern is a 
function of the angular coordinates 6 and m, while the incremental area dA 
is a function of p and c|), i.e. pdpd0. Thus, it is necessary to determine 0 
and CO in terms of p and <}), It can be easily shown that 

P / 

cos? + r sin? cos(<fi-iti) 

COS0 = » — (2a7) 

Vl + (P/h)^ 
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Figure 2A1. Geometry for the 'flat-sea impulse response computation. 
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and 


-tan^ + (-^) cos 
n 

COSO) = ' 

+ tan^C - 2(P/h) tan^ cos(4’-'t^) 


(2A8) 


It should be noted that the (f>-dependence always appears as the difference 
thus we can ignore <j> since the integration on (j> is over a complete 2tt 
radians. The first set of antenna patterns recorded by GE (ref. 2A4) indicated 
that the vertical port antenna pattern was nearly independent of co (circularly 
symmetric) out to about “20 dB. Revised antenna patterns (ref. 2A5) , however, 
show that both the horizontal and vertical ports exhibit asymmetric pattern 
shapes. When the antenna pattern is asymmetric, the m-integration in (2A6) 
is very difficult to perform and must be done numerically. As a result, 
one tends to lose the physics of the problem in mathematical detail. In 
order to simplify the problem, it was assumed that the antenna pattern was 
indepdendent of the azimuth angle coordinate co and of the form* 

_ sin^e 

f.( 0 _e ^ ( 2 A9 ) 


The constant y determines the beamwidth of the pattern and should be chosen 
to yield an average of the two principal plane beamwidths when dealing with 
asymmetric patterns. Substituting (2A9) in (2A6) and using (2A7) yields 



(1 + e ) 


2 2 2 2 1 
(cos K + Esin2^ cos<J> + e sin 4 cos (J>) / d(|) 


(2A10) 


*This approximation is assumed valid out to the -10 dB points in the 
pattern. 
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where e = p/h and has been replaced by since the integration is over 
2 tt radians. Denoting the <|>-integral as Q and rearranging terms yield 


Q = e 


1 - 


1 

cos ^ 

1 + 


+ b 


2tt 

/ 


acoscb T bsin <|) ,, 
e d(? 


where 


a = 


4e sin2^ 


(1 + e ) 


b = 


,2 .2, 

4e sin i 


(1 + E ) 


Since 


- bsin^0 


n=0 


. -.n , n . 2n 
(-1) b sin p 

( 


n! 


equation (2A11) can be integrated exactly and 


r 2 " 

4 ^ cos K . 1 . 

- 7 1 2“ > 

^ L (1 + e ) J / ■, 


Q = 2/tT e 


n=0 


(-1)" r(n + |) 
r(n + 1) 


(■^)" I 


where the I (•) are Bessel functions of the second kind and order 
n 

2b/a = Gtan^, Q becomes 




cos g 

(1 + E^) 


n=0 


(- 1 )" r(n + |) 
r(n + 1) 


• I 

n 



sin2^ \ 
(1 + E^) / 


(2A11) 


(a) 


n. Since 


(etan^)^ 


(2A12) 
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It should be noted that since both ^ and p/h are very small, the series in 
(2A12) will numerically converge very rapidly. Substituting (2A12) and 
3 = p2/h2 in (2A10) yields the following expression for the flat sea impulse 
response, 




vV r 6^t - /I + 6^ 


L (4tt) 
P 


^ f 

7t)\^ Jo 


W 


1 + 3 


exp 


i-+- 


cos^C „ . 2, 


00 

E 

n=0 




r(n + 1) 


d6 


(2A13) 


Using the property of the 6-function, the integral in (13A) is equal to the 
integrand evaluated at 


t - — /I + 3 = 0 or 3 = c^t^/4h^ - 1 

c 

and 




3 2 

(4tt) h L 

P 


a°(^o> 




COS K 

2 2 2 
c t /4h 



• 2. 
sin K 


OO 


■ E 

n=0 


(-1)" r(n + i) 
r(n + 1) 




1 


in 


-1 tan^ 


Converting to the two-way incremental ranging time, i.e., r = t - 2h/c, the 
flat-sea impulse response becomes 
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(-1)^ r(n + |) 
r(n + 1) 



(2A14) 


The total mean, return is a convolution of (2A5) with (2A14) , i. e. 


o 


C2A15) 


where the lower limit in (2A15) is taken to be 0 since P^(t,^) is identically 
zero for t < 0. In general, (2A15) must be integrated numerically; however, 
for C = 0° (nadir) some simplification is possible. In particular, the flat- 
sea impulse response is, for ^ = 0“, 


P^(t,0) 


r .2 ,2 

ttG X 
o 


3 2 

(4tt)Vl 




Since zzr ~ << 1 for nadir operation, 
2h 


itG^ X 

(t,0) w 

(4TT)\^Lp 





(2A16) 


and since the antenna beainwidth is so narrow we may take ~ 0° • Sub- 
stituting (2A16) and (2A5) in (2A15) yields the following expression for 
the mean return waveform on nadir. 
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7TG X 
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3 2 

2(4tt) h L 


ft)-- 
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(0) exp 


A(t-IOO)^ + “ 


2c 

yh 


- ACt-100) 





erf 



- A(x-lOO) 

7r 



(2A17) 


where 


A = 


1 - 


(l H- 


It should be noted that /a is proportional to the pulse length and not to the 
square root of the pulse length; this arises because of the way in which the 
pulse shape was defined (see equation 2A1) . 

The two principal-plane antenna patterns are shown in Figure 2A2 
for vertical port excitation. As can be seen, the Gaussian pattern 
approximation provides a reasonable fit*, however, the parameter y must 
be changed by about a factor of two between the m = 0“ and 90® pattern 
cuts in order to maintain the fit. Figure 2A3 is a plot of the mean 
return waveform as computed from (2A17) with 


K 


V 


^2 ,2 
ttG X 
o 


2(4tt) W 


The difference between the two curves for y = 7.88 x 10 and 3.66 x 10 

is 0.88 dB at the peaks of the two mean returns. These two curves provide 

bounds on the true mean return which will lie somewhere betx^een the y = 

—A “A 

7.88 X 10 and 3.66 x 10 curves. Also shown in Figure 2A3 is a mean 

-3 

return with y = 5 x 10 which is just about midway between the other two 

curves. Hence, it was assumed that the true vertical port antenna pattern 

could be represented by the Gaussian function given in (2A9) and with y = 
-3 

5 X 10 . Unfortunately, time was not available to reduce GE*s antenna 
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Figure 2A2. Vertical port principal plane antenna patterns (In Figure 2A1 
» = 0 corresponds to the y = 0 plane and a; = 90° corresponds ’ 
to the X - 0 plane) . 
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Figure 2A3. Nornialized nadir mean return for various antenna pattern tapers. 




pattern data and compute <Prs 3 (t. 0 )> these data; however, the circularly 

symmetric Gaussian pattern with y - 5 x 10 should be a goot. computational 
approximation. 

Figure 2A4 is a plot of the normalized mean return waveform at nadir for 
various values of rms waveheight These curves clearly show that both the 

shape and amplitude of the mean return waveform are dependent on onJ^ 
through variations in for reasonable values of i. e. ^ ^ meters.* In 

deriving equation (2A17) , we assumed that could be taken from inside the 
integral since the angle = tan"*^(p/h) would not differ appreciably from the 
pointing angle of the antenna. Using the dependency of on as determined 
in ref, Al, equation (2A15) was numerically integrated for various pointing 
angles of the antenna out to 15.6“. There was no appreciable difference between 
computed with cr°(<|^) inside the integral and <(P computed with 

1T3 S ITS S 

a taken outside of the integral. Figures 2A5 through 2A9 are plots of the 

mean return for antenna pointing angles of 0“ and 0.5“ , 3“, 8“ , and 15.6°, 
respectively. These results were computed by numerically integrating (2A15) 
using (2A14) and (2A5). 

Figure 2A10 shows how the peak-of-the-mean return decreases as the antenna 
is scanned off nadir, excluding effects. This curve shows the decrease in 
the peak-of-the-mean return as a function of pulse shape, scan angle and antenna 
pattern. The total decrease in the peak-of-the-mean return as a function of 
scan angle is obtained by adding the variation of (in dB) to this curve. For 
a wind speed of 5 knots, 0° changes by about 17 dB as the pointing angle varies 
from 0“ to 15.6°. Thus, the total change in the peak-of-the-mean return from 
0° to 15.6° is about 36 dB. This is the dynamic range which would be required 
of the receiver if it detected the peak-of-the-mean return. 


*When C = 2 meters, the significant waveheight, H is 8 meters. 

Tn ' 
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Figure 2A4. Normalized nadir mean return as a function of the rms waveheight ^ . 
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Figure 2A5. Normalized mean return for antenna pointing angles of 0® and 0.5® off nadir 
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Figure 2A6. Normalized mean return for an antenna pointing angle of 1.5® off nadir. 






ts 



T/^cReM6NTr‘^L ^ 7^0 'Way Dclav //m6 ^ 27 ( A//c/^os£coA/i?s) 


Figure 2A7. Normalized mean return for an antenna pointing angle of 3® off nadir. 




Figure 2A8» Normalized mean return for an antenna pointing angle of 8® off nadir. 
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APPENDIX 2B 


SIMULATION OF THE AGC r-FACTOR 

This section presen :s the results of our investigation of the relation- 
ship between the "mean of the peaks" and the "peak-of-the-mean" or the r- 
factor. The need for this relationship comes about because of the manner 
in which the Skylab S-193 altimeter operates. (The GEOS-C AGC system operates 
similarly in the acquisition mode.) This AGC scheme basically involves 
development of a control voltage which is a function of the peak signal 
voltage encountered during an AGC gate interval. Since data processing 
procedures depend on the peak of the average waveform instead, there is an 
obvious need for a connecting relationship. 

The first work relating to r-factor, by Lambert (ref. 2B1) , comprised 
a state variable simulation of the r-factor values and use of closed-form 
detection statistics from Rice (ref. 2B2) . The present effort was under- 
taken expecting to corroborate Lambert’s findings; the results, however, were 
found to be in substantial disagreement with his results. We argue the 
validity of our data based on the added degree to which the S-193 receiver 
is modeled. For example, rather than assuming Rayleigh or Chi-Squared 
detector statistics, the detector/video amplifier statistics are actually 
generated to the extent to which a valid system model exists for the S-193 
hardware. 

A block diagram of our hybrid computer simulation is shown in Figure 2B1. 
All computations prior to and including the video amplifier were performed 
using an EAI 380 analog/hybrid computer and the reraaining operations utilized 
a PDP8/e minicomputer. Receiver thermal noise was not included in this 
simulation and the results are felt to correspond to SNR values greater than 
10 dB. Other simulation details are identical to those given in our 1972 
report (ref. 2B3). Three 100 ns waveform cases were computed corresponding 
to 0, 1, and 1.5 degrees off-nadir; a few sample values of the mean wave- 
forms are plotted on the theoretical waveforms* in Figures 2B2 and 2B3 and 
the correspondence is seen to be quite good. Time did not permit computation 

*These waveforms are slightly different from those computed in Appendix 
2A. TTjey do, however, correspond to the waveforms employed in Lambert's 
computations. This was done in order to study the exact same situation which 
Lambert studied. 
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Figure 2B2. Comparison of theoretical and simulated mean return waveform at nadir. 






Figure 2B3 


Comparison of theoretical and simulated mean return waveform at 1.5® off nadir. 


of other, large off-nadir angles. 

Figure 2B4 shows a comparison of the two r— factor results. Whereas 
Lambert’s result decreases rapidly over the angle range of 0-1.5 degrees, 
our result displays an r— factor of 0.69, 0.58, and 0.53 at angles of 0, 

1.0, and 1.5 degrees. Note that the discrepancy between the two curves 
is considerably larger than the expected angular dependence of o . 

In order to verify our simulation results on the r— factor problem, 
we have investigated the AGC calibration data presented in ref. 2B4. Figure 
2B5 is a plot of the AGC (voltage) analog-to-digital converter output versus 
altimeter receiver input power as measured at the output of the antenna. 

These data were obtained from the flight model of the altimeter and were 
selected because they represent normal operating conditions, i.e. 25 C, 

10 MHz bandwidth and a 600 ns AGC gate window. Unfortunately, ref. 2B4 
did not state what the input power was; however, it must be the mean-of-the- 
peaks of the input power waveform because of the manner in which the AGC 
voltage is measured. 

The input pulse shape was triangular with a 100 ns rise time and a 500 ns 
decay time. This shape is a very good approximation to the mean return pulse 
shape at nadir (see Figure 2B2) , The "clean pulse" curve was obtained with 
"just the pulse input, and the input power in this case is equivalent to the 
peak-of-the-mean return power.* The "clutter pulse" curve was obtained with 
the input being a product of noise and the triangular pulse, which is approxi- 
mately what the actual sea-scattered signal will be at nadir. In this case, 
the input power will be equal to the mean-of-the-peaks of the input signal 
since the AGC is a "peak-picker" device. For the input power between -95.5 
to -60.5 dBm, we see that the clean pulse curve requires 3 dB more input power 
than the clutter pulse curve in order to yield the same AGC output. From 
our previous reasoning, we note that the mean-of-the-peaks of the return power 
is thus equal to the peak-of-the-mean return power plus 3 dB. Therefore, for 
an input power range of —95.5 to —60.5 dBm and at nadir, the r— factor is —3 dB. 
Comparing this result with Lambert's computations and the simulation results 
in Figure 2B4, we see that the ^3 dB agrees very well with the simulation result 
of 3.24 dB (0.69 in voltage ratio). We therefore conclude that the r~factor as 
derived from our simulation effort is substantiated by the data shown in 
Figure 2B5, of off-nadir angles of 3, 8, and 15.6 degrees is advisable. 

*Since the input waveform is deterministic, the peak of the mean detected 
signal is equal to the mean of the peaks which is what the AGC measures. 
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Figure 2B5. AGC calibration curve obtained by G.E. 
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Another interesting factor is given in Table 2BI below. 


Table 2BI. r-Factor Versus Number of V>faveforms Processed 

Numbe r r-factor 

1 
2 
4 
8 

16 
32 

These data correspond to an off-nadir angle of 1.5 degrees. The results 
show that an approximate steady— state condition should exist for the S-193 
AGC time constant (~50 samples). 

In summary, we feel that the comparisons shown identify a serious data 
processing information deficiency in regard to o° measurement with the 
S-193 system. To solve this problem, we recommend that an r-factor study 
be conducted using the S-193 breadboard hardware in conjunction with 
pulse shaping circuits which approximate computed waveshapes over the 
range 0, 1.5, 3,, 8, and 15 degrees off-nadir. 


.9992 

.8761 

.5892 

.5436 

.5588 

.5491 
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APPENDIX 2C 


COMPUTER PROGRAMMING FOR TRACKER AND SAMPLE-AND-HOLD SIMULATION 
Hardware Details for Waveform and Sampling Simulation 

The analog computer used to generate the simulated ocean-scattered 
radar pulses was a EAI-380, and the 12-bit minicomputer used to simulate 
either the altimeter sample-and-hold gates (for waveform experiments) or 
the altimeter ' s ■ splitgate tracker was a PDP8/e. The following list 
provides additional details of the hardware available for waveform and 
tracker simulation: 

Available Laboratory Hardware 

1 - PDP8/e Computer 

1 - Type KM8-E Memory Extension (total of 8K memory) 

1 - Type KE8-E Extended Arithmetic Element 

1 - Type ADOIA A/D Converter with AH04 Sample and Hold 

Amplifier 

2 - Type AA50 D/A Converter 

1 - Type DR8-EA 12-Channel Buffered I/O 
1 - Type DK8-E Real-Time Clock 

(Above two items available but not used in the 
simulation described here.) 

1 - EAI 380 Analog Computer 

Storage scope, X-Y plotter, etc. 

The relationship between the analog and the digital computers is 
presented in Figure 2C1. The analog computer is operated in a repetitive 
pulse mode of operation and the square-wave output from its time base is 
used as a synchronization or timing signal to one of the A/D channels 
of the digital computer. This sync signal will be denoted VS (the wltage 
for ^nchronizing) in the following discussion. 

After satisfying synchronization criteria, the digital computer samples 
the simulated radar return pulse signal, denoted SIC in the following. After 
sampling a preset number of return radar pulses, the results are displayed 
on an oscilloscope or are printed out on the computer's Teletype unit. 

The digital computer programming was done in the PDP8 assembly language 
to produce a large, general-purpose sampling subroutine which could handle 
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Figure 2Cl. Relationship between Analog and Digital Computers. 








either the altitude tracker or the waveform sampling functions as desired; 
this subroutine is called from FOCAL, the BASIC-like conversational 
language available on the PDP8. We used FOCAL to handle the details of 
formatting, I/O, and so forth, but used the assembly language sampling 
subroutine for its speed and efficiency and because of the specific control 
required for timing details. 

General Sampling Subroutine Structure 

To discuss our sampling subroutine, we use the terms gates , weights , 
and individual samples whose relationship is shown schematically in 
Figure 2C2. All timing is relative to the zero-going transition of the 
sync signal, VS, and the three important relative delay parameters dj^, 
d^, and are also indicated in the figure. The entire sampling sequence 
is started a variable delay time, d^, from the sync signal's zero tran- 
sition. There are a total of N gates, with each gate consisting of a 

G 

total of N individual samples. The individual samples are multiplied by 
w 

different weights, W, and the result for a given gate is a weighted sum, 
(sample- x W- + sample^ x w + sample x w ). a delay, d , between each 

X A A JM^ X 

weighted sample within a gate determines the total gate width. At the end 
of each gate a delay, d^, determines the between- gate spacing. 

The delays d^, a^, arid 'd;^ are d.mpTemented by a counting loop which, 
in effect, uses the basic computer cycle as a clock. The gates are 
all equal in width, equally spaced in time, and have identical weighting 
within the gate; this entire sequence is "free-running" when started at 
d^ after the sync zero. The basic computer timing cycle thus provides the 
relative spacing for the entire sampling sequence, with the sync details 
providing the absolute time reference. 

In our sampling subroutine, the maximum number of gates is 128, with 
each gate comprising from 1 to 128 individual samples - consequently there 
is a maximum of 128 weights. At the end of each sequence of _< 128 gates, 
the gate results update a set of 128 separate gate sum locations, so that 
the program accumulates a running mean. (As a minor point of interest, we 
also store a set of ^ 128 squares of individual gate results; the gate 
variances can be obtained, if desired, from the sum of squares and the 
square of the sura of individual results.) 
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Figure 2C2, Time Relationships in Sampling Subroutine. 



Uses of Che Sampling Subroutine 

To use the sampling subroutine for a simple waveform sample, we set 

N ^ 1, set weight =1, adjust the delay parameters d and d as desired, 
w 2 . ^ . 

and set N for anywhere from 1 to 128, to obtain a set of point samples 
G 

on the waveform as averaged over Ng sample cycle (that is, over Ng successive 
simulated returns) , 

For use as a split-gate tracker, we set = 2, set = 128, all weights 
to 1/128, adjust d^ as appropriate, set d^ so that the gate width is correct, 
and set d^ as appropriate for the time— separation of the two gates. In 
order to use the sampling subroutine for closed— loop tracker simulation, 
the two gate results must be transferred back to the main FOCAL program 
where the operation "first gate minus half the second gate" is performed, 
and the closed loop transfer characteristics are used to calculate the new 
position of the tracker. This position is set into the sampling subroutine 
as a new value of the delay-f rom-sync parameter d^^. 

Finally, the use of the sampling subroutine to simulate a set of eight 

sample-and-hold waveform sampling gates is straightforward; set = 8, 

adjust widths by d,, set d^ = 0 (to obtain immediately adjacent gate results), 

and set the N and the individual weight values as appropriate to the 
w 

characteristics of the sample-and-hold gates. If the sample-and-hold acts 
as a perfect integrator, all weights would be equal. If the sample-and- 
hold acts as a final-value-only hold, all weights would be zero except for 
the last which would be 1. A case intermediate between these extremes would 
be the sample-and-hold having a RC-type of behavior, and in this case the 
weighting would be maximum at the final weight and exponentially reduced 
from this at progressively earlier weights. 

Computer Program of Sampling Subroutine 

There is little point in reproducing program listings for the sampling 
subroutine just described. The listing is not only machine-dependent but 
also would be confusing since the subroutine is part of a general package 
which includes a floating point variable storage and recall area, an integer 
storage and recall area, and a routine to repetitively display on the oscillo- 
scope (via the D/A converters) the contents of the integer storage area. The 
interrelationships of these functions and their relationship to the overall 
FOCAL system plus standard 8K FOCAL overlay have led to several different parts 
of the sampling subroutine being located in various discrete positions in the 
memory of the computer. 
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While a complete listing is not provided, an overall fl.w chart for the 
sampling subroutine is shown in Figure 2C3. The delays are identified 
merely as block.s, whereas actually the delay operation involves counting up 
in a loop. Input parameters include the delays DL, DI , and DB 
d^ respectively) , information specifying channel and gain for the A/Ds for 
both the sync signal VS and the waveform channel SIG, the sync threshold 
VT, and the desired total number of complete sample cycles NS. Three 
separate subroutine entry points are indicated, and their differences are 
as follows: Entry 1 loads the complete set of parameters and the full set 

of sample weights, sets all results to zero and starts the sampling; Entry 
2 sets results to zero and starts the sampling for parameters and weights 
as already loaded; and Entry 3 permits acquiring additional sample cycles 
for a new higher value of NS. The actual number of complete sample cycles 
acquired is variable NT, and the results are stored in the array RESULT(I), 

I - 1 to NG, where NG is the number of gates. (Although not shown on the 
flow chart, an additional array RESULT2(I) stores the sum of squares of 
individual results for variance calculations.) Two loops are invloved in 
the details of synchronizing with VS; the first loop checks to see that VS 
exceeds the threshold VT, and when this is satisfied the second loop waits 
until VS falls below VT which occurs on the sync signal’s trailing edge 
transition to zero. 

One additional function indicated in the flow chart is provided by 
variables MAX and TMAX. TMAX stores the highest individual sample point 
obtained for each single sample cycle (during which a total of NG x NW 
individual sample points will have been acquired) , and MAX contains the 
sum of the TMAX values. This provides information needed in an investi- 
gation of the peak-of-the-mean vs, mean-of-the-peaks problem, the r-f actor 
problem described in Appendix 2A of this report. 

On subroutine exit, the results RESULT(I) [and RESULT2(I)], NT, and 
max are available to the main FOCAL routine. All subsequent processing 
and print-out is done there. 


Discussion 

The sampling subroutine as described has been implemented in a "real- 
time" situation. We indicate real-time in quotation marks, because the 
radar pulse times and bandwidths have been scaled to fit the capabilities 
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of the available analog and digital computers. One could consider situations 
in which the digital computer input might be magnetic tape containing already- 
digitized data. A straightforward FORTRAN program could be written to 
accomplish the general operations indicated in Figure 2C2. The flow chart of 
Figure 2C3 would be changed only at the details of synchronization, and by 
writing "READ SIG” in place of "Sample SIG (signal A/D)" and "READ VS" in 
place of "Sample VS..,". The delay operations would be different in detail; 
rather than counting in a loop, the "Delay, DL" operation would consist of 
skipping a number (related to DL) of digitized samples from the input tape 
(and similarly for delays d^ and d^) - 



CHAPTER 3 


TASK E, GEOS-C SATELLITE PROGRAM 
3.0 INTRODUCTION AND SUMMARY 

This chapter presents the results from investigations conducted under 
Task E, the contract task which is related to the GEOS-C satellite program. 

A summary of the contents of each subsequent section of Chapter 3 is provided 
immediately below. 

Altimeter data processing objectives and requirements are examined in 
Section 3.1 in terms of signal-processing characteristics of the presently 
planned altimeter hardware. Data considerations related to the acquisition 
of altitude data, geodetic slope data, and waveform information are discussed 
relative to basic characteristics of the altimeter data which derive from the 
on-board altitude tracker and the effective illuminated area. 

Section 3.2 gives estimates of the signal-to-noise ratio (SNR) under which 
we expect the equipment to operate, and presents recommendations for system- 
level tests relative to SNR. 

Section 3.3 considers telemetry data rate requirements necessary to support 
the presently planned experiments. We find that low and high data rates of 300 
and 90.00 bits /sec are needed.* 

The final section, 3.4, reviews presently planned calibration and qualifi- 
cation tests and gives recommended additions. 

An appendix to this chapter discusses the finite radar spot sizers smoothing 
effects and the consequent limitations on surface feature measurement; a possible 
computational program to predict such effects is described. 

3.1 Altimeter Data Considerations 

The General Electric Company (GE) , the GEOS-C hardware contractor, is 
presently planning to design the altitude tracker to have a closed loop band- 
width of approximately 4 Hz, for use in all modes except acquisition. The 
resulting filtering due to the tracker loop is to be followed by an on-board 
averaging process which will also affect statistical properties of the altimeter 
data. Autocorrelation properties of the tracker output and the tracker-plus- 
averager output are shown in Figure 3.1 for inputs that are uncorrelated. 


*At this time, it is planned to telemeter only 10 sample-and-hold (S&H) 
gate values. If all 16 S&H values are telemetered, the required data rate is 
approximately 14,000 bits /sec. 
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Figure 3.1. Autocorrelation functions for tracker output and averaged tracker output. 


(Derivation of these quantities will be given in 3.1.1.). Examination of 
Figure 3,1 shows that the correlation span is correspondingly longer for 
the averaged output,* and the detailed calculations given in 3.1.1 show 
the averaged output data will be essentially statistically independent 
over 0.1 second intervals. Also note that present plans are to telemeter 
the tracker error signal in synchronism with each received pulse. Thus, 
in principle, the altitude error signal for each pulse will be available 
for post-flight data reconstruction; this reconstructed signal should be 
free of most on-board tracker constraints. 

The selection of on-board processing characteristics should be subject 

to the following considerations: 

1) For altitude measurements only , there is a consensus of opinion that 
averaging altitude data points over a one— second time period will be neces 
sary to yield satisfactory variance properties in the geodetic data. On 
the basis of designing to a particular altitude noise level per unit time, 
it makes little difference whether the averaging is accomplished in the 
hardware or on the ground. This statement of course assumes that the on- 
board tracking and averaging devices are statistically near— optimum. 

2) A second (and generally misunderstood) requirement derives from the 
need for obtaining geodetic slope data. Note that a slope measurement is 
akin to a dif ferentation process, and the data processing techniques would 
be expected to differ from those of altitude measurement. The work of Cohen 
and Zondek (ref. 3.1) indicates that a considerably higher data rate may be 
needed for this purpose. It has been suggested that the maximum data rate 
of 100 altitude samples per second may in fact be best. This subject is 
somewhat complex; in a continuous — time analogy it represents estimation of 
the derivative of a stochastic process. In the discrete case, some forms of 
measurement noise can be shown to lead to a solution involving a minimum in 
the relationship between error and sample rate. This is in contrast to the 
usual practice of sampling a function at the highest practicable rate, a 
valid viewpoint for deterministic processes but not necessarily appropriate 
for stochastic processes. 

3 ) System parameters chosen to optimize data collection for either altitude 
or slope measurements, as in 1 and 2 above, are not necessarily best for 


*These results are somewhat different from those contained in the GE 
proposal, mainly because GE used a sin x approximation of tracker bandwidth 
effects. 
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purposes of waveform experimentations. We feel that a second, narrower 
bandwidth option should be available in addition to the presently planned 
4 Hz bandwidth. The reasons are: 

a) For waveform measurements this narrow bandwidth would be 
desirable to avoid having to attempt tracker motion reconstruc- 
tion with the attendant risk, telemetry problems and data 
processing complexity. With a 4 Hz bandwidth, tracker jitter 

will be approximately one half the pulse length (for moderately high 
SNRs); reduction of the bandwidth to 0.5. Hz, for instance, would reduce 
the variance by about a factor of three. This reduction should yield 
acceptable mean waveforms based upon on-board averages. 

b) The lower bandwidth value should constitute an aid to 
acquisition since, in the final step, the signal must be deter- 
mined to within ~20 nanoseconds (ns) from a 200 ns uncertainty. 

Detailed calculations for a Swerling II model show that from .5 
to 1 second of search time is required to produce acceptable 
detection and false alarm probabilities. Because of the altitude 
rates possible, the tracker should be able to maintain its slew 
rate (from the previous 200 ns tracking step) during the search 
interval. Therefore, a narrow tracker bandwidth would be 
desirable for acquisition. 

c) The lower bandwidth would permit acquisition of altimeter 
data in the global mode without further on-board averaging; it 
may also represent a better match to averaging characteristics 
of the footprint effect to be discussed in the next paragraph. 

Control loop calculations have been performed for a 0.5 Hz cross- 
over bandwidth based on velocity and acceleration values of 660 ns/ 

2 

sec and 6 ns/sec * respectively; the results show that the tracking 
error is ~1 ns. The corresponding velocity and acceleration 
coefficients are approximately .1000 and 10, These figures thus 
indicate that a 0.5 Hz bandwidth. Type 1 servo is capable of 
accommodating the target dynamics involved. 

4) In general, tracker bandwidth can be sho\m to be the determining 


*This acceleration is larger than the orbit derived value, and is 
based on differentiation of a curve fitted to the Puerto Rican Trench 
data of Von Arx. 
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factor in establishing short term correlation properties of the ensuing 
altitude data; however for studies in- which ocean-surface variations 
comparable to the illumin.-.ted spot size are of interest, tracker band- 
width effects are probably not the dominant process • Appendix 3A presents 
our preliminary work in this area; we argue that the tracker bandwidth 
determines the rate at which statistically independent estimates of the 
smoothed ocean surface process can be obtained, but that there is a 
smoothing function resulting from the (pulse-limited) spot size which 
may dominate the short term ocean surface autocorrelation characteristics. 
Quantitative values cannot be given at this time, and more work is definitely 
needed in this area. 

With these remarks on GEOS-C data acquisition philosophy, we turn to 
an analysis and discussion of altimeter instrument-determined data char- 
acteristics . 

3.1.1 Analysis of Altimeter-Data Variance Characteristics as Influenced 
by the On-Board Tracker .- This section presents results relating to the sta- 
tistical characteristics of the altimeter output data, under the assumption 
that each returned pulse represents a statistically independent quantity. 

The analysis parallels that given in the GE proposal (ref. 3.2) and page 
references from GE’s proposal are cited. The principal difference here is 
that the exact transfer characteristic of the tracker is used. We find that 
the presence of the on-board averager has an effect on the statistical 
properties of the intensive mode output data no different from that obtainable 
by a corresponding change in tracker bandwidth. 

With the closed loop tracker transfer function (given on GE’s page 

1. a. 2-160) as 


B(s) 


K (t s + 1) 

^1 

s^ + (1 + Kt_) s + 
1 
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using s = jo) where u is radian frequency, the magnitude squared function 
is 


= B(s) B*(s) 


2 2 2 
(Kt^) V + K 


t V + [1-2 KCt^-t^) + (Kt^)^] + K 


2 , 2 . .2 


In these equations, K is the loop’s velocity gain constant, denoted as K = K K 

* o 

by GE's pg. 1. a. 2-159. We next determine the autocorrelation by the 'transform 


P(t) = 2 ^ 


/ 


4>(w) cos(wt) dw 


which can be evaluated using the following integrals from page 411 of ref. 3.3 
for a>0, b>0, |t|<ir/2. 


Integral 3.733, No, 1, 


CO 

/ 


cos (ax) dx 


+ 2b^x^cos(2t) + b 


TT r -ab cos(t)l sin[t + 


ab sin(t)] 


(2t) 


Integral 3.733, No. 3. 


CQ 

/ 


X cos(ax)dx 


x^ + 2b^x^cos(2t) + b^ 


TT r -ab cos(t)l sin[t 

“ 2b J si 


>- ab sin(t)] 
in(2t) 


Using these, we get 

2 - ^ 1^/ 1 cos(t)J 


p(t) ^ 


0 ^ 
2/ K/t, 


sin(2t) 


sin[t - /K/t^T sin(t)] 


r— 

+ sin[t + /K/t-, T sin(t)] 

Kt i 
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where 


Since p (0) = 


1 -1 r ^ 

t - 2 2Kt^ 


(Kt^) 


2 1 


2/K/t~ sin(2t) 


the normalized autocorrelation function is 




- (cos t)rj ( sin[t - /K/t^ (sin t)T] 


P (t) = 
n 


^ 1^^2 

1 + — ^ sin t 


Kt, 


t /t 

+ •• sin[t + /K/t^ (sin t)T]] 


The function R^, a function of N, is defined as 


R, = 




1 N-1 ^ 

q=o k=o 


, where T = time between individual 
pulses . 


(See GE page 1. a. 2-60, where is defined as the "variance reduction factor 
due to the average of N tracking gate values to obtain average two-way range 
delay estimates.") To determine the effect of averaging on P^(t), K(t) for 
N pulses is defined by 


N-1 N-J. 


Kn(^) = “2 

N q=o 


E 


P^(| (T+q-k)T^| ) 


[clearly *^(0) = ^ 1 
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(This K,, is the same as GE's , pg. 1. a. 2-57, the "covariance function 

of the random error in the average two-way range delay estimate. ) 

These equations have been evaluated for a closed loop bandwidth 
of 4 Hz with the results as given in Table 3.1; graphed values were already 
given in Figure 3.1. Note that the tracker autocorrelation is <0.1 for 
T ^ 0.1 sec and that the averaged autocorrelation is .173 at x = .1 sec. 
Suppose the tracker bandwidth is reduced to approximate the averaged output. 
In the system variance term 



with . 1 

sec averaging 

= 1.71* 

without 

.1 sec averaging 

= 1.0 

with averaging 

= 4 Hz 

without 

averaging 

= 4 -f 1 


and to first order the variance is identical in the two cases. We there- 
fore conclude that the 0.1 sec averager is superfluous. 

For the one-second average of the intensive-data mode, the variance is 



AT 


R 



+ 



+ 


2 

""^Ql 


8.04 

16.53 


+ 


.22 

10 


+ 3.26 = 


(1.95 ns)^ 




or 29 cm standard deviation in height, where the values from GE page 1. a. 2-62 
and -68 were used and was recomputed for N = 100. This result shows the 
dominance of the last term; the standard deviation in height would be 27 cm 
if only this last term were present. 


*This compares with the GE value of •^1.78 obtained by approximating the auto- 
correlation function as a sin(x)/x. 


4 * 

'According to conversations with E. L. Hofmeister of GE, the 
be one- fourth the value given in the proposal. 



term will 
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TABLE 3,1 


TRACKER AUTOCORRELATION FOR NO AVERAGING 


AND .:'OR 0.1 SECOND AVERAGING 


Time x 

Autocorrelation 

- 

in seconds 

Tracker Output Only 

Averaged Tracker Output 

0.00 

1.000 

1.000 

.01 

0.851 

0.983 

.02 

0.711 

0.939 

.03 

0.583 

0,871 i 

.04 

0.469 

0,786 1 

.05 

0.366 

0.688 1 

.06 

0.276 

0.582 1 

.07 

0.197 

0.473 1 

.08 

0.131 

0.366 1 

1 

.09 

0.075 

0.264 

.10 

0.028 

0.174 1 

1 

.11 

”0.009 

0.098 ! 

.12 

-0.038 

0.035 j 

.13 

-0,060 

-0.015 1 

.14 

-0.076 

-0.055 

.15 

-0.087 

-0.085 ! 

.16 

-0.094 

-0.107 ! 

.17 

-0.097 

-0.121 1 

.18 

-0.097 

-0.130 1 

i 

.19 

-0.094 

-0.133 

.20 

-0.090 

-0.133 

.21 

-0.084 

-0.129 

.22 

-0.078 

-0.123 

.23 

-0.071 

-0.116 

.24 

-0.064 

-0.107 

.25 

-0.056 

-0.097 

.26 

-0.049 

-0,087 

.27 

-0.042 

-0.076 

.28 

-0.036 

-0.066 

.29 

-0.030 

-0.057 

.30 

-0.024 

-0.045 
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3.1,2 Discussion of Altimeter Variance as Influenced by Digital 

Delay Gate Start Signal Quantization Error .— Next, we turn to a discussion 

2 . • 
of the term o in GE’s tracker varxance equation 

^Q1 



The quantities here are: 
o 

= variance of the random error AT in the average two-way range 
estimate, 

2 

a = variance of the random zero-mean tracking error, 

T 

= variance reduction factor due to averaging N tracking 
gate delay values , 


N = number of successive tracking gate values averaged for the 
two-way range estimate, 


2 

a 



variance of digital delay gate start signal quantization 
error, 

variance of digital delay gate stop signal quantization 
error. 


Approximate numerical values for the tracker variance from the above 
equation and based upon GE's proposal (except for the values which are 
based on the previous discussion) are as follows: 


^ + 3.26 (ns)^, for 0.1 second average (N=10) , and 

Ml 1.71 10 

R 


^ ^ h + 3.25 (ns)^, for 1.0 second averages (N=100) . 

ATj^ 16 . 5 100 

2 2 

The dominant terra in the 1 second average is the 3.26 (ns) from a ; 

2 ^01 

GE does not divide this variance by 10 (unlike the a term) 

Q2 
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because they state that, "...the quantization error is a slowly varying 

error. It may remain constant for time intervals on the order of seconds. 
Hence, ...assume it is fully correlated over the averaging time of N pulses, 
whereas ^ rapidly .rarying error which can assume any value in the 

allowed range on a pulse- to— pulse basis and hence can be assumed completely 


uncorrelated pulse-to-pulse. 

An extremely important question is the rate at which varies . 

Depending on that rate, the term may be different from the indicated 


3.26 (nanoseconds) 
tizing time step T 


This is because a 


^ has been derived from the quan- 


'Q1 


Q1 


under the assumption that the probability density 


function p(e^-,) is uniform over the quantizing range as shown below: 
Q1 





■Ql 


Figure 3,2. Probability density function for 


the error 


Under this assumption of uniform density, the mean and variance of the 
quantizing error are immediately found to be 


Q1 


CO 

■/ 


T 

2 


■ ®Q1 


= rv- 

J Q1 


"qi ‘*"qi 


and 0 


Q1 


00 

/ 




Tqi^/12 . 
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The same form holds for . For the intensive mode, the GE documents 

^02 

indicate = 6.25 ns and - 1*58 ns, leading to 


and 


„2 .16^. 3.26 (ns)2 


2 = = 0.22 (ns)^ 




But in order that be representative for altitude data requires 

eqi 

that the data be acquired over a long time span compared to the typical 

period of variation of in order that the assumption of uniform 

over the interval valid. Let correlation time of 

and we will consider three different cases depending upon whether 

is much larger than, much smaller than, or of the same order of magnitude 

as t , where t is the altimeter observation period of interest, 
o o 


Case 1. 


» t 


Q1 


Suppose it were possible to profile exactly the same ocean area a 
number of separate times, each profile taken in a time t^. Suppose that 
all factors except are constant. Then one could compare the m separate 

profiles obtained, and upon repeating the experiment for m large enough 
that (mt ) » T , the data would appear as in the following sketch, 

.. V 

Figure 3. J. 
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two-way 
range time 
estimate 



Figure 3.3. Possible data for a number of profiles of the same 


ocean area, for t 


r, >> t . 

°Q1 ° 


The point is that t does not change during any single profile, that 

2 

a would be related to the absolute altimeter accuracy, and that for local 

Q1 

geoidal feature mapping (over a time span t or less) the relative accuracy 

° 2 2 

of the profile is not affected by . That is, a_ .defined by T^. /12 is 

Q1 Qi 

related to the statistics of an ensemble of separate altitude measurements, 
and for practical reasons the actual altimeter will measure only one or two 
members of this ensemble of altimeter profiles. The manner in which 
should be interpreted depends on whether one's interest is in absolute or 
relative accuracy. 


Case 2. 


r << t . 
°Q1 


In this case e - can be treated the same as 


That is, the a 


IS corre 


^ 2 ^ 9 ^ 

ctly given by T^, /12 and the contribution of to o xn 


reduced by a factor of N since is decorrelated in a 

short time compared to the observation time. The variance formula would 
then be 


T T 

_£2_ _Q1^ 

12N 12N 


. 0.22 . 3.26 , .2 . ^ . 

H — TTj” + ”17;” knsi tor 0.1 second averages. 


In this case poses no particular problem. 
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Case ‘S. 


T and t of the same order of magnitude. 

'^Ql 

In this case any observed profile could contain one or several dis- 
crete jumps by in the estimated range T^. This means discrete altitude 

jumps of (6.25 ns, x (.15 m/ns, c/2) ^ 0.9 meters, and it is uncertain 

what the best procedure would be for removing this contamination from 

the altimeter data; possibilities range from simple data post-editing to 
time-varying Wiener filter implementation. It is worth remarking that 
changes of 0.9 meters in 0.1 second time intervals are small compared with 
the altitude rate of 100 meters/second which the tracker should accommodate. 

Based on the discussion of the above three cases, it is clear that more 
detail is needed about statistical characteristics of Cases 1 or 2 

would be strongly preferable to Case 3 , the intermediate case . If the time 
characteristic of falls into Case 3, it would be desirable to stabilize 
e . better in order to achieve more nearly a Case 1 situation. 

K(± 

3.2 Open Questions Regarding Estimated Values of Altimeter SNR 

In this section we review some of the uncertainties that arise in the 
“ca'lrcu'larijo'H.'' "bf '^r’e'Gelve d -pow e r . S u ch calculations mu st be regarded as 
approximations and engineering judgements and therefore used with due regard 
for needed safety factors. The three areas of uncertainty to be discussed 
are; first, what we have denoted the "threshold effect"; then, signal-to- 
noise-ratio (SNR) estimates; and finally, sea state experimentation. 

3,2.1 The Threshold Effect .— The threshold effect in the context of 
GEOS-C represents a possible degradation in system accuracy for low SNR 
values, which the radar equation does not. predict. The data to be discussed 
indicate that this effect might cause the intensive mode to be inoperative 
for moderate to large waveheight conditions. 

The threshold problem is one of considerable theoretical complexity; 
we have been unable to locate published results that are directly applicable 
to the GEOS-C problem. The literature also used the terms global, ambiguity, 
and anomalous errors to describe the threshold effect. This effect repre- 
sents one of the research areas of spread-spectrum techniques, and the results 
available in the literature are adequate only to identify the effect as a 
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potential problem area. The most recent textbook freatment of the problem 
is by Sakrison (ref. 3.4). He states (pg- 62), 

*‘Note that calculating the probability that the global maximum 
of I (s ,W) I falls outside the central interval is equivalent to 
solvini a first passage time problem for a Gaussian process, a 
classical unsolved problem. The problem of how to calculate 
estimates of mean-sqiiare error in the threshold region, or even 
to calculate for what value of signal to noise ratio ... threshold 
effects start to occur can be regarded as one of the important open 
problems in communication theory at the present time. 

Vfhile the following paragraphs present the results of our brief study of 

the effect, we feel that the existing theoretical work is not sufficiently 

developed to justify the considerable effort which would be needed to 

analyze the specific GEOS-C system; instead we recommend that acquisition 

and tracking tests be made, as a function of estimated loop losses, with 

the GEOS-C pulse compression and altitude tracker hardware as early as 

possible in the progress. 

Following are the results of our brief survey of the threshold effect. 
The paper which brought the subject to our attention is by McAulay and 
Hofstetter (raf. 3.5). In this paper, the authors focus on the threshold 
behavior of a chirp system through an analysis based on the Barankin bound. 

They find that for a time-bandwidth product of 100 and a SNR of 2.5 dB, 
the range variance is 7 dB larger than that predicted by the Cramer-Rao 
bound (which should correspond closely to the radar equation result) . They 
conclude that a linear 40:1 pulse compression system should operate at a 
SNR greater than 10 dB if meaningful unbiased range measurements are to 
be obtained. Two other important papers on the subject are by Seidman 
(ref. 3.6) and Glave (ref. 3.7). Both of these papers cast serious doubt 
on the utility of the Barankin theory. They show it to be a very loose 
and parametrically unstable bound which, for the cases considered, produces 
a very optimistic result. An abbreviated version of the Seidman result is 
shown in Figure 3.4, which indicates that the Barankin bound predicts 
performance in excess of the theoretical maximum (the classical Shannon 
bound). Both Seidman and Sakrison recommend use of the Cramer-Rao bound 
for the high SNR case and the Ziv and Zakai bound for the low SNR case. 
Seidman argues that the threshold can be located to within + 2 dB with 
this theory, for a relatively wide class of waveforms; this paper shows 
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3.4. Noisy channel performance bounds for a time-bandwidth product 



the threshold point to be between 14-18 dB for a time-bandwidth product, 
of 80. 

The papers thus far mentioned consider only point, non-fluctuating 
targets. To our knox^ledge, the only reference which considers other 
target characteristics is Van Trees (ref. 3.8). In this work, he considers 
a point Rayleigh-fluctuating target and derives error equations (p. 305) 
based on an approximate analysis. Figure 3.5 shows computed values based 
on his equations, for conditions we feel are close to the GEOS-C parameters. 
The Van Trees results are interesting; in this case a threshold behavior 
is encountered which is a function of the number of pulses available. 

Figure 3.5 shows this result and indicates that, for very low SNR values, 
a tracker which effectively operates on 10 samples (for example, a 5.0 Hz 
bandwidth and a prf of 100) will encounter range errors which are degraded 
by a factor of 10 compared to linear estimates. This theory involves a 
number of approximations and the results should be considered only quali- 
tative; however, the tracker-bandwidth tradeoff is most interesting. 

We discussed possible threshold effects in the altimeter during a 
telephone conversation with Mr. R, J. McAulay, and his comments and 
suggestions are as follows: 

1. He definitely feels that the altimeter pulse compression 
system will exhibit threshold behavior and is not aware of 
any theoretical work relating to the extended, fluctuating 
target problem. 

2. He recommends that the ability of the system to ’’maintain 
lock” be investigated via the Ziv-Zakal approach (point target) 
to determine if the global uncertainty exceeds the tracker "pull- 
in” range. As a consultant, he recommends Lawrence P. Seidman 

of Systems Control, Inc., Palo Alto, California. 

3.2,2 SNR Estimates .— The SNR value given in GE's proposal is ~14 dB 
for the intensive mode, based on a calm-sea value for o° of 9 dB. We 
estimate other system losses to be: 
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Figure 3.5 » Tracker rms error as a function of SNR and number of pulses 
averaged, fron Van Trees. 


antenna pointing error (1 degree off nadir) 

pulse compression sidelobe filter loss 

miscellaneous degradation 

(power output, noic e figure, atmospheric loss) 

Total 

Resulting SNR for calm seas 

Resulting SNR for SWH* 3-4 ft 

Resulting SNR for SWH <^6 ft 

Because of the marginal values of SNR shown, especially in 
to waveform data acquisition and to acquisition and maintenance 

lock, we recommend: 

1) The pulse compression sidelobe filter should be omitted. For the 
above SNR values time sidelobes would not be important, and also the system 
bandpass properties will provide some degree of sidelobe modification. 

2) The satellite altitude should be reduced by approximately 20% 
to achieve higher effective SNR. 

3) Again, system tests including pulse compression and the altitude 
tracker should be conducted as early as possible to document system 
performance for SNRs near 0 dB, particularly in the area of non-linear 
detection characteristics, 

■ 3 2.3 Sea State Experimentation .- Sea state measurement will degrade 
with SNR primarily because: 

1) The increased noise baseline will cause the signal component to 

be a smaller percentage of the recorded v/aveform. Therefore the uncertainty 
in measuring signal properties will be greater. 

2) At low SNR values the tracker variance will increase and S&H 
position uncertainty will become greater. Note that the S&H circuits will 
be tied to the range tracker and that time jitter (tracking error) will be 
superimposed on the on-board average waveforms. 

Ip, view of the potential degradation in range error due to the threshold 
effect. Category 2 is considered to be the most important error component, 
particularly since the S&H jitter is related to the high intensity mode 
tracker-bandwidth. For example, if the one second tracker jitter is .50 meter, 


*Significant Wave Height 


3.4 dB 
2.2 dB 

2 dB 

7.6 dB 

6.4 dB 

3 dB 
0 dB 

regard 
of tracker 
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the 10/second sample-rate jitter would be 1.5 meters or about one pulse 
length. These figures indicate that the considerable advantage of accumu- 
lating on-board average waveforms may be lost. Also note that previous 
studies have recommended that sea state experimentation emphasize rough 
sea conditions (see ref. 3.9). Present SNR estimates indicate this 
would not be advisable, even disregarding the threshold effect. The task 
of demonstrating the sea-state measurement concept is more difficult for 
small wave-heights. 

A problem associated with waveform studies is that of measuring 
signal fluctuation statistics. Recent data from the Wallops/NRL aircraft 
program suggest that present statistical scattering models may be of 
limited validity. The requirement for relatively clean individual wave- 
forms is more stringent here than in any of the experiments considered. 
Study may in fact show that this measurement is of limited usefulness 
for the SNR constraints encountered. 


3.3 Telemetry Requirements 

Tables 3,2 and 3.3 show the values which we consider to be minimum 
"ail'd — maxirnt 2 m“teieTnetry ‘ T-equirements for the various signal categories. 

These tables are largely based on the instrument data characteristics 
shown in Table 3.4 (which is taken from the GE proposal). The low data 
rate result shows that ~300 bits /sec are required unless the plateau gate 
is sampled every pulse. The plateau gate probably need not be sampled 
at this rate since the assumption of statistical independence of each 
sample and ocean surface homogeneity will allow a slower data accumulation 
rate (^lO/sec) . 

In the high data rate case, the' instantaneous plateau gate data are 
omitted since they are duplicated by the S&H data. Table 3.3 shows a bit 
rate of approximately 9000 bits/sec is needed if 10 S&H values are telemetered 
and a bit rate of • 14,000 bits/sec is needed if all S&H values are used. Note 
that the tracker error signal is the second largest item in Table 3.3. This 
arises since tracker reconstruction considerations dictate that this quantity 
be densely sampled. 
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TABLE 3.2 

LOW r vTA RATE TELEMETRY SUMMARY 


Bit Rate (BPS) 
Min Max 


I. Low data rate parameters 
Housekeeping quantities 
Approx. 6 values, 8 bits each 


Clock Data 


System Outputs 
same as I— mode 


(20 bits /minute) 
2/sec 


5 48 

(1/flO sec]) (1/sec) 

nil 

10 192 

(1/[10 sec]) (2/sec) 


II. Medium data rate parameters 

Status Words 8 bits, 1/sec. 

AGC 8 bits 

Altitude word upper 16 bits, 1/sec. 
Tracker error signal 8 bits,l/sec. 

III. High data rate parameters 

Altitude word 16 bits, 10/sec. 
Inst, plateau gate 8 bits 


8 8 

16 80 

(2/sec) (10/sec) 


16 16 

8 8 


160 160 

80 800 

(10/sec) (100/sec) 


Total 


298 1312 
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TABLE 3.3 


II. 


III. 


HIGH DATA RATE TELEMETRY SUMMARY 


Bit Rate (BPS) 



Min 

Max 

Low data rate parameters 



Housekeeping quantities (temp. & voltage) 



Approx. 6 values, 8 bits each 

5 48 

(1/[10 sec]) (l/sec) 

Clock Data (20 bits/minute) 

System outputs* - typical quantities: 


nil 

avg, noise power, avg. specular gate, 
avg. ramp gate, avg. plateau gate, calih. 
data, average waveforms, loop lock signal 
assume 12 parameters @8 bits 

10 

(1/[10 sec 

192 

]) (2/sec) 

Medium data rate parameters 



Status words 8 bits, l/sec 

8 

8 

AGC 8 bits 

16 

(2/sec 

16 

) (10/sec 

Altitude Word upper 16 bits, l/sec 

16 

16 

High data rate parameters 



Altitude Word 16 bits, 10/sec 

160 

160 

Tracker error signal 8 bits, lOO/sec 

800 

800 

Subtotal 

1010 

1304 

Waveform sample and hold data 100/ 

8000 

12,800 


(10 S&H) 

(16 S&H) 

Total 

9010 

14,104 


*Note that the instantaneous plateau gate is omitted in the I mode. 



TABLE 3 . 4 

GEOS-C ALTIMETER DATA OUTPUTS 
(From GE Proposal, Figure l»a.l-5) 


5MH/ REFERENCE 




IDO Hz SYNC 

CLOCK SIGNAL GNO A 
CLOCK SIGNAL GND B 


ALTITUDE OUTPUT COMMAND 

ALTITUDE SHIFT GATE 

ALTITUDE SHIFT CLOCK 

ALTITUDE OUTPUT 

ALTITUDE SIGNAL GND A 

ALTITUDE SIGNAL GNO B 


D1 COMMAND INPUT 

ECl COMMAND INPUT 

CK1 COMMAND INPUT 

COMMAND SIGNAL GND lA 
COMMAND SIGNAL BND IB 

D2 COMMAND INPUT 

EC2 COMMAND INPUT 

CK2 COMMAND INPUT 

COMMAND SIGNAL GND 2A 
COMMAND SIGNAL GND 2B 


SW1 STATUS VAORO -««* 

SW2 STATUS WORD 

SW3 STATUS WORD — * 

SW4 STATUS WORD --J 

SW5 STATUS WORD 

SW6 STATUS WORD 

STATUS WORD SIGNAL GND A 
STATUS WORD SIGNAL GND B 


TT1 THERMISTOR TRANSMITTER 

TT2 THERMISTOR RECEIVER — * 

TT3 THERMISTOR G-TRACKER 

TT4 THERMISTOR l-TRACKEH 

TT5 THERMISTOR WAVEFORM SAMPLERS 
TT6 THERMISTOR BIT/CAL-^ 


THERMISTOR BUS 


VlVf^) RELAY DRIVE VOLTAGE — 
V<AF)L0 FREQUENCE OFF SET — 
VR^l RECEIVER MIXER CURRENT 
V(RMI RECEIVER MODE 


V(Pj) TRANSMIT POWER 


'AGC 


RECEIVE SIGNAL LEVEL 


viNi ave. video noise POV/ER -« 

VIA/S) AVE.PVYR. IN ATTITUDE/SPECULAR GATE 

VqIRIAVE. RAR3P GATE ENERGY-G mode 

Vq(PJ ave. PLATEAU GATE ENERGY G MODE 
Vg(Pi) INSTANT. PLAT. GATE ENERGY-G MODE -« 
VgJTj) TRACKING LOOP JITTER G MODE 

V,(R) AVE. RAMP GATE ENERGY ! MODE 

VjlPfAVE. PLATEAU GATE ENERGY-1 MODE 
V,(Pi) INSTANT. PLAT. GATE ENERGY ! MODE 
V,|Tjl TRACKING LOOP JITTER I MODE 


ANALOG SIGNAL GND AA 
ANALOG SIGNAL GND AB 

signal GND A 

SIGNAL GND B 






+14 VDC BUS A 
+14 VDC BUS B 
POWER GND A 
POWER GNO B 


> A9 


> A6 


> A9 


GEOS-C 

RADAR 

ALTIMETER 


>-A9 


A13 

ARl 

A7 

A12 

A13 

A1 I 


ARl, A7, All. A12, A13, A18 

PSl 

A4 

A3 

A3 


> A7 


AI2 


PSl, A2. A3, A4, A7, A12 




PSl 
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3.4 Recommended Calibration Tests 

Altimeter system level electrical tests have been examined with 
reference to the test categories utilized in recent GEOS-C meetings. 

These categories are listed and defined below: 

A. Functional Tests 

These tests constitute verification of the proper sequencing 
and appropriate responses of all mode commands which can be exercised 
after the satellite is in orbit. As such, these constitute the most 
isolated or ’’buttoned— up” form of electrical test. Tests will be under 
the control of the TAMS system and altimeter input will consist of clean 
signals applied via waveguide or RF inputs. 

B. Subsystem Performance Tests 

These tests consist of data acquisition of the subsystem (e.g., 
tracking loop) and major component level (e.g., antenna), and of tests 
based on use of built-in test /calibrate functions. From an engineering 
standpoint, these tests will provide the most detailed information available 
on the system. 

C. System Performance Tests 

These tests constitute verification of the performance of the 
major system functions such as acquisition, altitude measurement, average 
waveform measurement, signal fluctuation probability density measurement, 
and drift-rate/bias tests on these parameters. 

D. Engineering Data/Calibration Data Tests 

These tests will culminate in a final calibration data book 
which will provide the prime data base for experiment investigators . 

The data base will represent results of tests defined in this category 
and all appropriate test data obtained from other test categories. 

In the efforts reported herein, we have attempted to place ourselves 
in the position of a GEOS-C investigator and to examine present test plans 
in terms of experiment requirements. On this basis our results apply only 
to categories B and D of the intensive mode. In presenting these findings, 
the information covered by E. L. Hofmeister at the 8 November 1972 Status 
Review meeting is given in Tables 3. 5-3. 9; our recommended test additions 
are given in Tables 3.10 through 3.12. These items are considered to be 
self-explanatory in view of previous discussions, and will not be Individually 
examined here. 
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TABLE 3.5 


SUBSYSTEM PERFORMANCE TESTS - GLOBAL MODE 


1. Transmitter 

Peak power 

RF f requency/pulse spectrum 
Pulse shape - single and burst 
PRF 

Burst pulse spacing 
Transmit pulse time jitter 
Transmit loss 

2. Antenna 

Performance measured at lower level ; antenna assembly . 

3. Receiver 

Noise figure 
Receive loss 

AGC/Tracking sensitivity 
Receiver dynamic range 
AGC dynamic range 
AGC dynamic range 
AGC linearity 

Receiver bandwidth; bandpass characteristic 
Detector; performance measured at lower level 
Video bandwidth; frequency response to AGC loop 
and waveform samplers 

4. Signal processor 

a . AGC loop 

Loop bandwidth; frequency response (acquisition and 
Loop step response tracking) 

b. Tracking loop 

Gate configuration; width, spacing, and spacing for burst 

Loop bandwidth; frequency response 

Loop step response 

Loop response to ramp input 

Track gate step size 

c. Altitude output 
Averaging process 
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TABLE 3.6 


SUBSYSTEM PERFORMANCE TESTS - INTENSIVE MODE 

1. Transmitter 

Peak power 
RF center frequency 
Pulse spectrum 
Pulse shape 
PRF 

Transmit pulse time jitter 
Transmitter phase shift 
Transmit loss 

2. Antenna 

Performance measured at lower level; antenna 
assembly 

3. Receiver 

Noise figure 
Receive loss 

AGC /Tracking sensitivity 
Receiver dynamic range 
AGC dynamic range 
AGC linearity 

'"'■'tercel ver "trandwidlih; bandpass characteristic 
Signal-to-noise ratio improvement 
Time sidelobe levels 
Compressed pulse shape 

Detector; performance measured at lower level 
Video bandwidth; frequency response to AGC loop, 
tracking loop, and waveform samplers 

4. Signal Processor 

a. AGC loop 

Loop bapdwidth; frequency response (acquisition 
Loop step response ■ and tracking) 

b. Tracking Loop 

Gate configuration; positions, widths, and spacing 

Tracking gate step size 

Loop bandwidth; frequency response 

Loop step response 

Loop response to ramp input 

c. Altitude processing 
Averaging process 

d. Waveform sampling system 
Sample gate width 
Number of sampling gates 
Sample gate spacing 

Reference time for sampling gates 
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TABLE '3.7 

SUBSYSTEM PERFORMANCE TESTS - GLOBAL & INTENSIVE MODES BUILT-IN TEST/ 

CALIBRATE 


1. Time delay calibration 

Instrument delay bias measurement; internal vs. 
external measurement 

2. Transmit/receive power calibration 

3. Receiver /AGC calibration 

Shape and level of IF signal //I 
Shape and level of IF signal //2 

4. Signal processor calibration 

Shape of video test waveform //I 
Shape of video test waveform if 2 


/ 


3-27 



TABLE 3.8 


SYSTEM PERFORMANCE TESTS 


Global Mode Tests 

GM 1 Acquisition Test 

GM 2 Height Measurement Test 

GM 3 Instrument Drift Rate Test 

GM 4 Instrument Bias (Calibrate/Test) Test 

Intensive Mode Tests 

IM 1 Acquisition Test 

IM 2 Height Measurement Test 

IM 3 Instrument Drift Rate Test 

IM 4 Instrument Bias (Calibrate/Test) Test 

IM 5 Average Impulse Response Test 

IM 6 Echo Amplitude Probability Density Test 
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TABLE 3.9 


ENGINEERING DATA/CALIBRATION DATA TESTS 

A. Final calibration data package consists of data from these tests 
plus selected data fx\)m performance tests 

B. Some important engineering/calibration data 


1. Transmit pulse shapes : photos of RF and detected RF 

transmit pulses 

2. Point target return pulse shapes: couple transmit pulses 

into receiver and photograph pulses at pre-compression 
IF test point, IF test point, and video test point 

3. AGC calibration curves: AGC control voltage versus peak 

input pulse power for different signal types and 
temperatures 

4. S & H amplitude calibration curves: S & H output voltage 

versus input voltage for each gate at several temperatures 

5. Calibration curves for other analog output signals: TBD 

6. All characteristics of built-in test/calibrate system 


• Calibrated attenuation value 

• Shape and level of IF test signals 

• Shape of video test waveforms 

• AGC voltage, tracking point, waveform sampler 
values ■•f ox opexat ion-on- test- -signals 

7. Probability density functions and autocorrelation functions 
for selected analog output signals for clean/noisy test 
signal inputs 

V 

AGC 

^(N) 


'^(A/S) 

Vg<R) 

Vg(P) 

Vg(P1) 

Vg(T1) 

lAWl 


V^(R) 

Vj(P) 

V^(Pi) 

V^(Ti) 

IIWl 


IAW16 


IIW16 
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TABLE 3.10 


RECOMMENDED ADDITIONS TO SUBSYSTEM PERFORMANCE TESTS* 

1. Transmitter 

a. Statistics of start of altitude countdown relative to transmit 
pulse (Eqj^) 

b. Transmit frequency stability 

c. TWT time delay versus H, V. level over regulation range 

2. Antenna 

a. Test to verify multipath free range conditions (e.g., demonstrate 
1/R signal dependence) 

3. Receiver 

a. Total receiver linearity test (r-f through video amplifier) 

b. Noise figure test of a batch of mixer/diodes to sample aging 
characteristics over several months and test environmental 
dependencies 

c. Measure receiver dynamic range with ACC voltage fixed (open 
loop) 

d. Measure pulse compression SNR improvement for post compression 
IF values of SNR in the range of 0-13 dB 

e. Measure receiver p re- and post-detection bandpass characteristics 

4. Signal Processor 

a. Measure tracker jitter statistics using a clutter signal with 

documented fluctuation statistics and for SNR (post- compression IF) 
values from "break-lock" to + 10 dB 

5. Subsystem performance tests/built-in test and calibrate 

a. Instrument delay bias measurement versus temperature, warm-up 
period, and accumulated instrument running time. 


*Most parameters should be measured as a function of expected variations 
in subsystem supply voltage, temperature, input signal level, etc. 
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TABLE 3.11 


RECOMMENDED ADDITIONS TO SYSTEM TESTS 


Performance Tests 

a. Noise figure with all subsystems in operation 

b. Acquisition time versus r-f signal level 

c. Tracker break-lock mean value 

d. Internal delay stability 

e. Transmit signal leakage (with r-f output separately 
terminated) relative to level and time delay of 
internal calibrate signals 
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TABLE 3.12 

RECOMMENDED ADDITION'^ TO ENGINEERING DATA/CALIBRATION DATA TESTS 


A. Measure pulse compression time sidelobes with reasonable time delay 
between transmit pulse to ensure that receiver recovery effects are 
not present (e.g., use a coaxial delay line or return from a corner 
reflector) . 

B. Simultaneously record the time difference between the on-board 
filtered altitude output and a sample-and-bold (S&H) position, 
with the tracker operating with a clutter signal, to determine 
whether or not S&H position reconstruction is necessary. 
Simultaneous records of the tracker error signal are also 
desirable. 


C. Measure drift characteristics of the tracking gate positions with 
respect to the transmit pulse. 

D. Measure receiver response to a short pulse (^one nanosecond) . 
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APPENDIX 3A 


As a number of authors have noted, the determination of geodetic data 
by spacecraft experiments is similar to the study of a noisy signal after 
passage through an electrical filter. For example, the geoid can be deter- 
mined from orbit perturbations of the spacecraft; in this case the equivalent 
filter is a low pass filter as the higher harmonics are not sensed by the 
spacecraft. The particular appeal of the radar altimeter is that it permits 
the direct measurement of higher order surface harmonics, at least over the 
ocean and to the extent that effects of ocean currents and meteorological 
phenomena are either absent or can be corrected. Even for this case, there 
is an effective higher frequency (shorter surface wavelength) cut-off which 
is related to the size of the radar footprint (among other things) . 

Typically, one calculates the radar footprint radius R (for pulse-length 
limited geometries) from the relationship 

R = (3A1) 

where h is the satellite altitude, c is the speed of light, and T is the 
radar pulse length; it is then often asserted that the footprint diameter 
(2R) represents the minimum surface wavelength which the altimeter can 
determine, and that the altimeter response is essentially flat for surface 
wavelengths greater than 2R and zero for wavelengths less than 2R. We feel 
that it is correct to approach this subject with the point of view. of deter- 
mining an effective transfer function to describe the satellite altimeter's 

response to the amplitude a of each surface wavenumber k (k = 2ir/X ) but 

s s s s 

that it is not adequate just to assume this transfer function flat out to a 

cut-off wavenumber k = 2ir/(2R) and zero for k > k . Instead it should be 

c sc 

possible to calculate the transfer function, and this appendix sketches the 
first steps in such a calculation. We have been limited by time and money 
but the problem is an important one which merits further study. 

The point of view used is as sketched below: 
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Note that the process labelled "Antenna filtering effects" implicitly 
Includes details of the transmitter and receiver impulse response, of the 
transmitted pulse shape, and of the altitude tracker processing. The 
starting point is the ocean surface’s description by a vector wavenumber 
spectrum, wavenumber Any given k could equally well be described by 

a magnitude k, a phase relative to the chosen (arbitrary) origin of the 
surface coordinate system, and a polar angle describing the orientation of 
k relative to the surface coordinate system's x~axis. We make the assumption 
that any given wavenumber does not change appreciably in direction over an 
ocean area centered at nadir and somewhat larger than the radar footprint . 
The ocean surface coordinate system’s x-axis is arbitrarily chosen to be in 
the direction of any specific wavenumber magnitude, k^ , being considered. 
Figure 3Al summarizes the problem's geometry and the notation to be employed. 

As can be seen from Figure 3A1, a non-zero ocean surface elevation , 

due to the wavenumber k , will cause a given incremental scattering area 

s 

located at p,<{) to be at a different distance, r, from the antenna than would 
be the case for a perfectly flat sea. Consequently the time at which this 
incremental area contributes to the radar return will be different from 
the flat-sea case, and the summation of return signal from the entire ocean 
-sur-f aGe:?iwouaxi..be £xpe.n^ to have a time-behavior different from the flat-sea 


value. 

This summation is made easier by noting that « h for any case of 
practical interest; thus to negligible error, the angle 6 is unchanged from 
rhe flat— sea value. Also, and co are the same as for the flat sea. The 

angle i{^, defined in Figure 3A1 as the angle between the z— axis and the straight 
line drawn from the satellite to the scattering area, will be taken as the 
angle that the radar pulse is incident upon the scattering area. This last 
assumption may be slightly risky as it ignores the slope of the local sea 
surface relative to the plane z = 0, but for the k^ values of interest to 
geodetic satellite programs the assumption is probably acceptable because of the 
small slopes anticipated. 

From Equation 2A6 of Appendix 2A of this report, the ocean surface 


impulse response is given by 


s (t,0 

a 

s 


Lp(4ir)^ 


j" <S(t - 2r/c) f^(0,m) 

scattering 
surface area 


dA 


(3A2) 
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z 



Assume that, over a subsatellite area 
somewhat larger than the radar "footprint 
the surface wave number 1^ is kg in magni- 
tude, is in the x-direction, and. has a 
phase angle as relative to x = 0. The 
amplitude of kg is ag, so the surface 
elevation at any point is given by Zg, 

z = a sin(k pcos(|) - a ) 
s s s s 


p = polar radial location of 
scattering area. 

ca - projection on z = 0 plane of 

the antenna pattern’s azimuthal 
angle . 

k = 27 t/A . 

s s 

h - satellite altitude 

r = distance from scattering area 
to satellite. 


i 

Figure 3A1. Satellite-ocean surface geometry, and summary of notation. 
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(For convenience of notation in later paragraphs, we use S rather than the P 
of Appendix 2A,) The various quantities not already defined in the present 
section are as defined in Appendix 2A, Notice that Equation 3A2 contains two 
approximations. The first is the use of the angle w to represent the antenna 
azimuthal angle whereas w is actually the projection of the azimuthal angle 
onto the plane z = 0. This is unimportant, particularly as one is so often 
forced to assume antenna patterns of azimuthal symmetry as was the case in the 
calculations later in Appendix 2A. The second approximation is the already- 
discussed use of ip in the backscattering cross-section , The integration 

over scattering area surface is carried out using the p,(f> coordinates, 

dA = pdpd0 , (3A3) 

and so we need a relationship between r and p for the 6-function in (3A2) . 

From the triangle whose sides are r, p, and (h-z ), find that 

s 

r h Vl + (p/h)^ - 2z /h , (3A4) 

s 


to an approximation good for all times of interest to the satellite problem. 
But r = ct/2 as well, and using the ranging time t defined by, 

t = + T , (3A5) 

we can find the relationship between p and t to be specified by, 

:?: ( |1 )^ + chT + 2hz . (3A6) 

^ s 

Notice that unlike the flat-sea case, the relationship between p and x is not 

unique because of the cj> dependence contained in the non-zero z term of (3A6) . 

4 ^ 

The r in the denominator of the integral 3A1 may safely be approximated by 

r^ = h^(l + p2/h2)^ . (3A7) 

Then assuming the antenna pattern f(0,tii) to be independent of lo and denoting 
the pattern by f(0), the integral (3A2) becomes 
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2tt 


'a = 2^ / \ 

s J s 

o 


, for T > 0 


(3A8) 


where 


S^(t,5,*) = 

o 

s 


2 2 
G X 
o 


2 4 

2(4tt) L h 

P o 


/ 


“ f^(6) a°(ip) 6[p“V(^)^ + chx + 2hz ] 

/ s 


pdp 


(1 + p^/h^)^ 


From Appendix 2A, we have 


. cos ^ + p/h sin C cos(c|)-(f)) 
cos o = ^ 


VT7 


(3A9) 


p /h 


From Figure 3A1, 


tan ip = p/h (3A10) 

and by earlier definition, 

z = a sin(k pcostf) - a ) (3A11) 

s s s s * 

Equation (3A8) above is written in terms of S (t,C »<{>); this notation is 
introduced here for later reference. While (3A8) as written above implies 
p-integration inside the <J)~integral , the order of indicated double integration 
could obviously be reversed and (3A8) rewritten to show (Jj-integration inside 
a p-integral. 

Thus (3A8) gives us the desired impulse response function for any specified 

k and with the p and <j> dependence of quantities in (3A8) specified by (3A9)- 
s 

(3All), but because of the complexity of the geometry and particularly because 
of the form of the 6-function argument for non-zero a^ in (3A11) , we cannot 
proceed any further toward the (desired) closed-form result. Expression (3A8) 
merely serves as a kind of symbolic shorthand to describe the two-dimensional 
numerical integration which has to be performed. Presumably a grid or cell 
size will have to be chosen on the surface with the 6-function providing 

a yes-no decision on whether each cell’s contribution is to be included in 
the running subtotal which is the digital computer's approximation to the 
double integral. A judicious choice will have to be made for the "width" to 
allow for the 6-function; because of the finite grid cell size and finite 
computer precision, the argument of the 6-function cannot be required to be 
exactly zero but must be provided a small width (either side of zero) within 
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which the 6— function will be non-zero. It is difficult to specify ahead of 

time how to choose these various parameters for a numerical approximation to 
(3A8) ; probably it will be necessary to vary the parameters to empirically 

ascertain where the numerical integration’s result will become reasonably 
independent of surface grid sizes and 6~widths. This procedure may be 
expected to be quite expensive in computer running times. 

Thus, while the result (3A8) may appear straightforward symbolically, 
it is of little use for practical reasons. We will now specialize the 
problem to ^ = 0, again assuming azimuthal symmetry of the antenna pattern. 
Figure 3A2 below reproduces just the plane r, p, h from Figure 3A1, at some 
instant of time x, for x > 0. - 


z 



Figure 3A2. Geometry for time x > 0. 

It is seen that r is the distance from satellite to scatterer which we 

would have if a =0, and n is the distance from satellite to scatterer for 
s 

non-zero a^. We calculate the time shift x^ - x corresponding to the differ 
ence between r^ and r by the following steps: 

r - r = - z cost|y (3A12) 

1 s 

cosi/j = ■ — . (3A13) 

Vl + (p/h)^ 

Since p = /chx , and z^ = a^ sin(k^/chx cos(|> - a^) , and since (x^^ - x) = 

2(r^ - t)/c, then for defined as ~ we find 
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(3A14) 


2a sin(a - k /cht cos(})) 
s s s 

A = — 

T C I p 

\1 + (ci/h) ' 

Here A is the time shift such that the impulse response function contribution 

T 

which would have occurred at t for a =0 now occurs at a new time t. , 

s 

T. - T + A . (3A15) 

1 T 

Now consider the use made of the surface impulse response function 

s 

For a transmitted signal P^(t), the received signal will be given by 

the following convolution 


(t,5) 

r t 


P,(x.O 





P^T-xp dx^ 


(3A16) 


If the properties of are such that 


P^“(f) = 0 


T < 0 


P (t) = 0 
t 


T > T 


then together with the property S (t) = 0, x < 0, the convolution expression 


becomes 


2it t 

d(f) 


f J Pt<x-Tp dx^ , X > 0 


and 


= 0 


x-T 


T < 0 


(3A17) 


In obtaining (3A17) , we substituted (3A8) for S and reversed the order of the 

X and (}) integrations. Notice that S (x,^,<t)) can be characterized (for C small 

1 

or zero) as describing two separate effects: 1) an amplitude effect due to the 

4 4 

change in geometry (the difference in ij/ and the 1/r vs l/r^^ as in Figure 3A2) ; 
and 2) a time shift effect as indicated in the discussion of the A^ defined by 
(3A14). 

The time shift effect will be very much more important than the amplitude 
effect for the case ^ = 0 and for the relative magnitudes of h and a^ in the 
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satellite problem. Notice from (3A8) that, if " G and K - 0> then 

S (t, 0,« is identical to S (t, 0). Thus for C = 0 and the already-assumed 
o o 

azimuthal antenna beam symmetry, a very good approximation to (3A17) is 
obtained by substituting for S^(t, 0,<J)) in the integrand, 

271 T+A^ 

Pr(T.O) f f So(Tj^.O) P^(t + - Ti) . (3A18) 

•b t+A -T 

T 

In (3A18) one can for all times of practical interest approximate (3A14) by 


2a 

A = — ^ sin(a - k /chx cos(^)) . (3A19) 

X c s s 

This expression for P (x,0) still will require fairly extensive computer times 

r 

for actual evaluation, but no unclear decisions or "widths of delta functions" 

y^qxiired as compared to our discussion of (3A8) • The virtue of (3A18) 

is that S (x,0) can be evaluated once and entered in the computer as a table 
o 

of values vs. x; it may even be possible to approximate S^(x,0) by a single 
function as in Appendix 2A, equation (2A16) of this report. 

If the problem is not restricted to C = 0, we see from (3A18) and (3A8)- 
(3A11) that for non-zero C, the 8-dependence on <j> will prevent the use of a 
previously calculated or tabluated set of values for (x ,^) in (3A18) . 
Instead,- the integral of (3A8) must be put into (3A18) , leading to a triple 
integration and the already-discussed problem of the proper choice of a 
"6-width". 

To see what (3A18) means let k 0 but keep a non-zero. Then 

s s 

2a 

A = — ~ sina , (3A20) 

T c s 


k ^ 0 
s 


and 


Pj.(t,0) = 

k -> 0 
s 


x+A^ 

J S(t^,0) P^(x + A^ - x^) dXj^ 


x+A -T 

X 


(3A21) 
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By changing variables, (3A21) may be rewritten as 


2a 

p (x - sina ) - 

r c s 

k ^ 0 


/ 

r-T 


SoC^ijO) P(t - dx^ . 


(3A22) 


s . . 

Thus for 0, the result is a simple time-shift of the signal from 

the flat-sea situation, with the maximum time shift being 2a^/c and the actual 
time shift determined by the phase angle which, in effect, described the 
surface harmonic phase at the satellite nadir point. The condition k^ ^ 0 
actually means k << 27r/(2R) with R as given by (3A1) . (Notice also that 
(3A22) becomes the usual time convolution ifa=0.) 

' Holding a fixed at some typical value (perhaps a couple of meters) , it 
s 

will be necessary to consider the details of the altitude tracker loop and 

determine the balance point for the different shaped waveforms P^(x,0) produced 

as k is increased. The entire process can be carried out for a number of 
s 

different k values, holding a fixed and setting a = tt/2 (the maximum 
s s s 

deviation from flat-sea) , and the results will eventually lead to a graph of 


the approximate form shown below. 
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This curve is the desired result, the w«%venuniber transfer function -which 
describes the effects of the footprint filtering process on the wavenumber 
spectrum. Of course, this discussion has been in the idealized limit (only 
waveform ensemble averages have been considered) , and so the results of a 
program of calculations based on (3A18) would represent the fundamental limit 
of the altimeter measurement of surface harmonics. We have not carried out 
such a program, but have described in this appendix how to begin on the 
problem. 



